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Abstract 

There is a common description of different intrinsic geometric flows in two dimensions 
using Toda field equations associated to continual Lie algebras that incorporate the de- 
formation variable t into their system. The Ricci flow admits zero curvature formulation 
in terms of an infinite dimensional algebra with Cartan operator d/dt. Likewise, the 
Calabi flow arises as Toda field equation associated to a supercontinual algebra with odd 
Cartan operator 8/86 — 68/ dt. Thus, taking the square root of the Cartan operator 
allows to connect the two distinct classes of geometric deformations of second and fourth 
order, respectively. The algebra is also used to construct formal solutions of the Calabi 
flow in terms of free fields by Backlund transformations, as for the Ricci flow. Some 
applications of the present framework to the general class of Robinson- Trautman metrics 
that describe spherical gravitational radiation in vacuum in four space-time dimensions 
are also discussed. Further iteration of the algorithm allows to construct an infinite 
hierarchy of higher order geometric flows, which are integrable in two dimensions and 
they admit immediate generalization to Kahler manifolds in all dimensions. These flows 
provide examples of more general deformations introduced by Calabi that preserve the 
Kahler class and minimize the quadratic curvature functional for extremal metrics. 
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1 Introduction 



The subject of geometric flows witnessed rapid development in recent years. It only fair 
to say that many areas of mathematics underwent rapid development as new tools be- 
came available for addressing long standing open problems - predominantly in differential 
geometry - via appropriately chosen evolution equations, also known as geometric flows. 
At the same time, several cases of such flows arose independently in physics with appli- 
cations ranging from classical mechanics to general relativity and quantum field theory. 
Roughly speaking, there are two general classes of deformation equations in differential 
geometry, namely intrinsic and extrinsic curvature flows. The former refer to geometric 
evolutions of Riemannian metrics g on a given manifold M that are driven by the Ricci 
curvature, in various forms, whereas the latter refer to geometric evolutions of subman- 
ifolds embedded in M that are driven by their extrinsic curvature. As such, they all 
correspond to dynamical systems in superspace, which is the infinite dimensional space of 
all possible metrics on a given manifold. The geometric flows provide systems of parabolic 
equations, being first order in the deformation variable t that is typically called time, 
and they are second or higher order in the space variables. However, the equations are 
quite difficult to solve in all generality, due to non-linearities. Although the short time 
existence of solutions for a given initial metric is guaranteed by the parabolic nature of 
the equations, their convergence to canonical metrics after sufficiently long time has only 
been analyzed under various conditions, and in connection to the problem of formation 
of singularities along the flows. 

All geometric flows share some common qualitative features with the linear heat 
flow equation, which has the tendency to dissipate any temperature perturbations all 
over space after sufficiently long time. Thus, curvature perturbations around canonical 
metrics tend to wash away after infinitely long time, provided that no singularities are 
formed along the way and that the corresponding geometric deformations are taken in 
some normalized form, so that the volume of the space remains invariant in time. It is 
precisely for this reason that continuous flows to constant curvature metrics have been 
particularly successful for exploring the relations among the geometric and topological 
structures of manifolds and prove various geometrization conjectures in low dimensions, 
as in the classic uniformization problem of Riemann surfaces. However, there is no 
general dictionary at this moment that connects geometric evolution equations to physical 
problems, apart from some special examples that are under intense investigation in recent 
years. It is quite plausible that field theory, classical as well as quantum, may provide 
a unifying framework for all different species of geometric flows by giving appropriate 
physical interpretation to the deformation variable t and to the structure of their driving 
terms. 

The Ricci flow is the prime example of this kind as it describes the beta function 
equations of two dimensional sigma models to lowest order in perturbation theory. The 
Calabi flow, on the other hand, makes its appearance in general relativity while studying 
the general class of four-dimensional Robinson- Trautman radiative metrics. Both these 
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flows correspond to intrinsic geometric deformations and they will occupy most of the 
present work. In fact, it proves advantageous to treat both of them in parallel, as we 
uncover some formal relations among the two, which exist in two dimensions and general- 
ize quite naturally to Kahler manifolds in all dimensions. Thus, some known techniques 
for the integration of Ricci flows will be found to admit immediate generalization to 
the Calabi flow, using the appropriate mathematical framework that will be developed 
in due course. As byproduct, the formal construction of solutions describing spherical 
gravitational waves in vacuum will be achieved in the context of four- dimensional gen- 
eral relativity, thus providing the main physical application of this work. Although our 
analysis is only applicable, as it stands, to intrinsic geometric flows, it is worth noting in 
passing that other important classes of evolution equations, such as those correspond to 
extrinsic curvature flows, might become tractable by similar methods. Such generaliza- 
tions are mostly interesting for a variety of problems in boundary quantum field theory, 
but they will not be included here. 

The primary aim is to develop some new algebraic techniques for casting intrinsic 
geometric flows in two dimensions into zero curvature form and prove their integrability, 
in a certain sense, using infinite dimensional Lie algebras. These results should be con- 
sidered complementary to the traditional techniques used in the mathematics literature 
for studying general properties and qualitative features of geometric evolution equations. 
Although strictly limited to two dimensions, where most of the algebraic techniques are 
becoming available, the results might also be of more general value given the embedding 
of two-dimensional manifolds into higher dimensional spaces that deform geometrically 
in one way or another. The present description views the geometric flows as Toda field 
equations associated to infinite dimensional algebras that incorporate the deformation 
variable t into the system of their defining commutation relations. This approach works 
well for the two-dimensional Ricci flow, but it also generalizes quite naturally to the Cal- 
abi flow using a formal relation between the two via super-evolution. In fact, as will be 
seen later, one follows from the other by taking the square root of the time evolution op- 
erator, and, thus, their formal integration by group theoretical methods appears to have 
many similarities. This relation will be further generalized to induce an infinite hierarchy 
of higher order evolution equations and associated infinite dimensional Lie algebras that 
enable to cast them in zero curvature form. 

The algebras used for the purposes of the present work deserve further study as they 
fall in the general (but poorly understood) class of infinite dimensional Lie algebras with 
infinite growth. Two-dimensional integrable systems based on infinite dimensional Lie al- 
gebras are quite interesting to study on their own, in all generality, and in comparison to 
the better understood cases of integrable systems that admit zero curvature formulation 
with finite dimensional gauge algebras. A key point in our analysis is the uneven treat- 
ment of the space coordinates and the time t, which is fully encoded into the structure 
of the infinite dimensional algebra used for the effective two-dimensional description of 
the equations. Thus, the systems appear to be integrable in the two-dimensional sense, 
whereas at the same time exhibit dissipative behavior in time without reaching a contra- 
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diction of terms. Finally, the algebraic structures underlying these flows might be of more 
general value if could be extended to encompass other interesting dynamical phenomena 
in geometrical theories with relevance to the physical world, such as gravitation. 

Actually, the framework that is developed here, based on Toda field equations for 
infinite dimensional algebras, admits a smooth limit when t-independent solutions of the 
two-dimensional Calabi flow are investigated. This class of solutions is also of paramount 
importance in general relativity as they describe Robinson- Trautman space-times of a 
certain algebraic type. In this will be seen in detail, there is a zero curvature 

formulation of the equation using Kac's K 2 simple Lie algebra with infinite growth. The 
results that we derive will also be compared to other works on this subject and they 
will be further used to construct formal power series solutions in terms of free fields as 
in ordinary Toda systems. As byproduct, a new description of the K 2 algebra will be 
obtained within the general framework of supercontinual Lie algebras that are introduced 
here for the first time and used for the zero curvature formulation of the Calabi flow, its 
variants, and descendants. Thus, the systematic study and classification of such novel 
algebraic structures seems to be a valuable problem of common interest in physics and 
mathematics. 

Some generalizations to geometric flows on Kahler manifolds of arbitrary dimension 
will also be considered. In that general context, the hierarchy of geometric flows provide 
special examples of a more general deformation problem posed by Calabi for minimizing 
the quadratic curvature functional, while preserving the Kahler class of the metrics. The 
standard Calabi flow is the most elementary example in this class of deformations with 
known physical applications. All of them, however, are naturally associated to Calabi's 
variational problem in the space of all possible metrics on a given Kahler manifold with 
Euler-Lagrange equations satisfied by extremal metrics, such as metrics with constant 
curvature. Thus, they have equal importance in geometry. 

The structure of this paper is organized as follows. In section 2, after briefly re- 
viewing the general set up of the main equations, the emphasis is subsequently placed 
on deformations of two-dimensional geometries, where the new algebraic techniques are 
becoming available for their (formal) integration by group theoretical methods. In sec- 
tion 3, we give account of the notion of continual Lie algebras, their main properties 
and the systems of Toda field equations associated to them in zero curvature form. The 
general method of integration is also illustrated for general choices of Cartan operators 
and further generalizations to the class of supercontinual Lie algebras are introduced. 
In section 4, applications of the algebraic method to the two-dimensional Ricci flow are 
briefly summarized together with the formal construction of the general solution in terms 
of free fields by Backlund transformation. In section 5, the method is extended to the 
two-dimensional Calabi flow by introducing an infinite dimensional algebra whose struc- 
ture is dictated by the super-evolution operator that connects it to the Ricci flow. In 
this context, it is more natural to view the zero curvature formulation of the Calabi flow 
as providing two-dimensional integrable system for a supercontinual Lie algebra, where 
the generators depend on the coordinates (t, 9) of i? 1 ' 1 superspace. This reformulation 
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also proves useful for the description of the general solution in terms of free fields, as 
for the Ricci flow. In section 6, the relation of the two-dimensional Calabi flow to the 
theory of spherical gravitational waves in vacuum is spelled out in detail using the class 
of Robinson- Trautman metrics that represent gravitational radiation in the exterior of 
bounded sources in four space-time dimensions. Then, genuine solutions of the Calabi 
flow correspond to type II gravitational backgrounds in Petrov's classification of space- 
time geometries. Type III metrics, which correspond to time independent solutions of the 
Calabi flow, are studied separately in section 7 by Toda field theory techniques and the 
connections to Kac's K 2 algebra are spelled out in great detail. In section 8, an infinite 
hierarchy of higher order geometric flows is constructed in two dimensions by iterating the 
formal relation between the Ricci and Calabi flows, via super-evolutions, and their zero 
curvature formulation is studied in all generality. In section 9, this hierarchy is shown to 
admit immediate generalization to Kahler manifolds in all dimensions and the resulting 
evolution equations provide special examples of a more general problem posed by Cal- 
abi for deforming arbitrary Kahler metrics toward extremal metrics. Finally, section 10 
contains our conclusions and a brief outline of some directions for further work. 

2 Intrinsic geometric evolutions 

In this section we present an account of the Ricci and Calabi flows, as they are defined 
in Kahler geometry, and then reveal a formal relation between the two that plays impor- 
tant role later. Although our results are mostly confined to one complex (or two real) 
dimensions, a more general, but brief, introduction is necessary in order to understand 
which of them can be extended to arbitrary dimensions. 

2.1 Preliminaries 

Let M denote a complex n-dimensional manifold, which for all practical purposes is taken 
to be connected and compact without boundary. Assume also that M admits a Kahler 
metric g, which is locally expressible in the form 



using a system of holomorphic coordinates z a and their complex conjugates z a with 
a = 1,2, • • • ,n. The Kahler condition on the metric is expressed by the existence of 
a locally defined, real valued, smooth function K(z,z), which is not unique, such that 
g a l(z,z) = d a d b K(z ) z). We also consider the associated real valued exterior form of 
degree (1, 1), 



which is closed and determines the principal cohomology class of the metric, [uj] . As such, 
it is an intrinsic invariant of the complex analytic structure on M. We are interested 



ds 2 = 2g al dz a dz b , 



(2.1) 




(2.2) 
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only in the space of Kahler metric in a fixed cohomology class; other Kahler metrics of 
the same class can be obtained by the transformation 

d 2 u 

9ab = 9ab + fcaQzb ' 

where u is a globally defined real valued function on M. 

For later use it is convenient to recall the expression for the Ricci curvature tensor, 
in complex notation, 

_ g^ogjdetg) 

and the Ricci scalar curvature, which is 

R = = g al R ai . (2.5) 

Note that the normalization used here is one half of the usual value found in most texts. 
Finally, the volume element is 

dV(g) = J n] = {^l) n det(g al )dz l A dz 1 A ■ • • A dz n A dz n (2.6) 

and it is also an invariant of the complex analytic structure on M. 

The Laplace-Beltrami operator for a Kahler metric g on M is defined as usual, using 
the analyst's sign, 

= \g-VN 3 <P = 9 al V a %<p = <f*-^ . (2.7) 

The integral J Atp dV(g) over the entire space vanishes by Stoke's theorem, since it has 
been assumed that M is compact without boundaries. Also, for later use, it is convenient 
to recall the definition of the fourth order operator, D L , which is an elliptic differential 
operator that acts on scalar functions as follows, 

D m a™ 8 (n b5 9 (n 9 (a 1 ™ ^ ^ 



dz a V dz b y in dz cKi> dz m > j J 

= AAn^% + V%. (2.8) 

It is a strongly elliptic operator that coincides with its complex conjugate Dl if and 
only if R is constant and, furthermore, satisfies the relation D L R = D L R for any Kahler 
metric. 

The geometric flows that will be considered next provide deformations of the Kahler 
metrics within a given cohomology class [u]. As such they can be formulated as parabolic 
evolution equations for the single function u(z,z;t) that arises in the class of transfor- 
mations (2.3), or equivalently Co = uo + ddu. However, we will follow the standard 
presentation based on evolution equations for the components of the metric itself, which 
are driven by the intrinsic curvature of M in various forms, depending on the type of 
flow, and which characterize the order of the non-linear differential equations. 
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In all cases, the metric is prescribed at some initial time, say t — 0, and evolves 
accordingly to later times. Short time solutions always exist, due to the parabolic nature 
of the equations, but their long time existence is not guaranteed and depends on various 
conditions and on the dimensionality of the manifold M. When they exist, however, have 
the tendency to converge to canonical metrics defined in appropriate ways. The advantage 
of using Kahler geometry is that all information about the curvature is encoded into a 
single function, the scalar curvature R, as in two dimensions. However, many important 
details also depend on the dimensionality of the Kahler manifold and its properties. 

2.2 Ricci flow 

The Ricci flow is the simplest example of a geometric deformation, which is driven by 
the Ricci curvature tensor associated to a metric g on the manifold M and assumes the 
following form (see, for instance, [1, 2], and references therein) 

dtQal = -Rab ■ (2-9) 

Actually, this equation can be defined on any Riemannian manifold of arbitrary dimen- 
sion, but we only consider here its application to Kahler geometry. More general, the 
Ricci flow takes the form 

d t g^ = -R^u + V M &, + V„^ , (2.10) 

thus also taking into account the effect of arbitrary (time dependent) reparametrizations 
associated to a vector field £. It is a parabolic second order non-linear differential equation 
that arose independently in physics and mathematics for different reasons. 

In physics it has definite meaning within the renormalization group analysis of two- 
dimensional sigma models, as it provides the rate of change of their target space metric - 
viewed as generalized coupling - with respect to the logarithm of the world-sheet length 
scale, [3, 4]. The Ricci curvature tensor arises as the beta function of the target space met- 
ric, to lowest order in perturbation theory of the quantum theory, and drives the evolution 
of the generalized coupling from the ultra-violet to the infra-red regime. Solutions to this 
equation do not correspond to conformal field theories, but they rather describe transi- 
tions between different closed string vacua in the weak gravitational regime of the theory 
induced by tachyon condensation. In this context, reparametrizations in target space that 
are associated to gradient vector fields, £ M = V M 0, can be interpreted as having a dilaton 
field in target space coupled to the metric. Likewise, fixed points of the generalized Ricci 
flow equation (2.10), satisfying the Ricci soliton equation = 2V M V;,0, correspond to 
non-trivial conformal field theories. The Euclidean black-hole in two-dimensional target 
space with the shape of a semi-infinite long cigar provides the simplest example of such 
a Ricci-Kahler soliton with axial symmetry. 

Another interesting application arises from the possibility to embed the Ricci flow 
into Einstein equations using two extra dimensions with coordinates u and v, as follows, 

ds 2 = -2dudv + g^ix; u)dx^dx u , (2.11) 
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provided that the light-cone coordinate u is identified with t and there is also a dilaton 
field linear in v coupled to the higher dimensional metric, [5]. This prescription allows 
to construct gravitational backgrounds - though of special type - for consistent string 
propagation, using solutions of the Ricci flow equation in two dimensions lower. Then, 
the renormalization group flow becomes the profile of a gravitational wave in space-time, 
and the geometry interpolates between fixed points of the flow in transverse space as one 
traverses the wave. In effect, this construction makes the Ricci flow to appear as special 
case of Einstein equations, with the appropriate ansatz, and it is the closest one can get 
to a similar derivation of the (two-dimensional) Calabi flow from Einstein equations, as 
will be seen later. 

In mathematics, it is often appropriate to use the normalized Ricci flow, which is 
defined as follows, [1, 2], 

< R > 

d t g^ = -R»„ + ^—^9^ , (2.12) 
where < R > is the average (mean) scalar curvature 

Im r [9) dv (g) , ,a\ 

S M dv{g) 

which is independent of the coordinates on the manifold and can only depend on t. The 
compensating term is designed so that the volume of the space M is preserved under 
continuous deformations, unlike the standard (unnormalized) Ricci flow that changes it in 
time. Thus, the normalized flow is more useful as tool for proving various uniformization 
theorems in geometry, since its solutions do not become extinct after some finite time, 
but they have better chance to exist for sufficiently long time and converge to canonical 
(constant curvature) metrics that correspond to the fixed points of equation (2.12) on M. 
However, the two variants of the flow are related to each other by performing rescaling 
of the metric and time reparametrization according to the rule 

g^ = f(t)g^, t = Jdtf(t), (2.14) 

which sets < R > (t) = (dimM)f'(t)/f 2 (t). 

In Kahler geometry, in particular, the Ricci flow has been employed for studying the 
existence of Kahler-Einstein metric on manifolds under appropriate technical conditions. 
It has already been applied successfully to two dimensional Riemann surfaces of genus 0, 
1 or g > 2 by re-deriving the classic uniformization theorem of Poincare that asserts the 
existence of constant curvature metrics with values +1, or —1, respectively; as a result, 
these surfaces are quotients of S 2 , R 2 or H 2 by a discrete subgroup V acting freely and 
isometrically. In higher dimensions the situation is promising but more complex. 
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2.3 Calabi flow 



The Calabi flow, on the other hand, which is only defined for Kahler manifolds, assumes 
the following form, [6, 7], 

9t9al = ^ ■ (2 - 15) 
It is volume preserving deformation without need for adjustment, unlike the Ricci flow, 
and it also preserves the Kahler class of the metric. The Calabi flow is a parabolic 
equation for the components of the metric, or equivalently for the Kahler potential, 
but it is fourth order in the variables z and z and, thus, more complicated to analyze in 
general. Critical points of the flow are called extremal metrics and they clearly encompass 
constant curvature metrics, if they exist on a given Kahler manifold. In this respect, the 
Calabi flow is used as a tool for studying the conditions for Einstein-Kahler metrics in 
geometry, and in conjunction with their possible obstructions. 

There is a variational problem in the space of Kahler metrics that allows to formulate 
the problem in general terms. In particular, consider the quadratic curvature functional 

S(g) = [ R 2 [g] dV{g) (2.16) 

J M 

defined for all possible metrics on M with fixed cohomology class [lu] . By computing the 
variation of S(g) with respect to the flow (2.15) one finds 

d t S(g) = f (-2RAAR - 2RR ai R ai + R 2 AR) dV(g) , (2.17) 

JM ^ ' 

since the associated scalar curvature flows as 

d t R = -AAR - R' a ~ b R ah . (2.18) 

Next, it is useful to employ the fourth order operator D L to cast the result of the variation 
into the form 

d t S(g)= [ (-2RD L R + 2RR' a R ja + R 2 AR) dV(g) (2.19) 

JM ^ ' 

and observe that (R 2 R' a )^ a = 2RR' a R a + R 2 AR. Thus, dropping the total derivative 
term, it follows immediately that 

d t S(g) = -2 f RD L R dV(g) . (2.20) 

JM 

The operator Dl is positive semi-definite and self-adjoint operator acting on functions 
on M. It can be further decomposed as Dl = L*L, where L is a second order operator 
mapping complex valued functions to holomorphic tangent vector valued (0, l)-forms, 

L<p = ip^g cl -^®dz b , (2.21) 
and L* denotes its adjoint. Then, using for simplicity the Hermitian product 

(<p, V)= / <p(z,z)$(z,z)W(g) , (2.22) 

JM 
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we have in total 



d t S(g) = d t (R, R) 



2(R, D L R) 



2(LR, LR) < , 



(2.23) 



showing that S(g) is a Lyapunov functional that decreases monotonically along the Calabi 
flow. 

The critical points correspond to extremal metrics on M, [6, 7], which according to 
the previous discussion satisfy the Euler-Lagrange equation DlR = 0, or equivalently 
the simpler equation LR = that describes its global solutions in a compact manifold M 
without boundary. This, in turn implies that the extremal metrics, when they exist, are 
characterized by the condition that their curvature is such that R a d a is a holomorphic 
vector field. In general, very little is known about the space of extremal Kahler metrics, 
but it is expected to have unique extremal metric in each Kahler class up to holomorphic 
transformations, apart from some special examples. Further analysis of the problem 
shows that if there exists a constant scalar curvature metric on M, then it achieves 
the absolute minimum value of S(g). In this case, all extremal metrics have constant 
curvature and the absolute minimum value of the quadratic curvature functional that 
follows by application of Schwarz's inequality, 



is attained. However, there can be situations where extremal metrics exist on M but 
they do not have constant scalar curvature; in fact, constant scalar curvature metrics and 
extremal metrics with non-constant scalar curvature do not co-exist in a single Kahler 
class. 

Whenever the existence of a critical Kahler metric can be guaranteed, then the so 
called Futaki-Calabi obstructions determine the necessary and sufficient conditions for 
the existence of a constant scalar curvature metric in Kahler geometry, [8, 7]. Last, 
but not least, there are counter-examples of Kahler manifolds that do not admit any 
extremal metrics and, hence, S(g) has no critical points; for, if extremal metrics exist 
they will necessarily be symmetric under a maximal compact subgroup of the holomorphic 
transformation group, but, in some cases, this fails to exist, [9]. Given these facts, it 
appears that the situation is rather complex in general and we will not attempt to fill in 
the technical details as they are beyond the scope of the present work; for an overview, 
see, for instance, [10, 11]. Here, we will be mostly concerned with two-dimensional spaces 
with the topology of sphere, which admit a constant curvature metric - that of the round 
sphere - and permit the solutions of the Calabi flow to exist for sufficiently long time 
and converge to the canonical metric, [12, 13]. The conditions for global existence and 
convergence of the Calabi flow to constant curvature metrics is a hard problem, in general, 
in higher dimensional spaces. 

The physical relevance of the Calabi flow is not yet known for manifolds with ar- 
bitrary number of dimensions apart from the lowest dimensional case of one complex 
(or two real) dimensions. It is quite remarkable that the two-dimensional Calabi flow 



S(g)> 



(J M R[g] dV(g)f 
J M dV(g) 



(2.24) 
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characterizes completely the solutions of four dimensional Einstein equations that de- 
scribe gravitational radiation from bounded sources, using the general class of Robinson- 
Trautman metrics for space-times that admit a geodesic, shear and twist free diverging 
null congruence. In this context, t is identified with the retarded time coordinate in 
space-time, as will be explained in detail in section 6. 



2.4 A formal relation between the two flows 

There is a curious relation between the Ricci and Calabi flows on Kahler manifolds of 
arbitrary dimension that manifests by squaring the time evolution operator. In particular, 
taking the time derivative of the Ricci flow, we have 

0p9ab = -Qj. R al = d a d b -^(log detg) , (2.25) 

which in turn can be evaluated by noting <9f(log detg) = g ab d t g a i = —g ab R a l — ~R- 
Thus, in Kahler geometry, deformations of the metric flowing a la Ricci always imply the 
equivalent relation 

gp9ab = -d a d b R . (2.26) 

Clearly, if the second derivative of the metric with respect to the Ricci time is identified 
with minus its first derivative with respect to the Calabi time, 

d 2 d . . 

w R = ~w c • < 2 - 27 > 

the two flows will be formally the same. 

The procedure above looks similar in vain to the way of extracting the square root of 
the Schrodinger equation in supersymmetric quantum mechanics for a free particle (for 
a recent account see, for instance, [14], and references therein). In the latter case, this is 
practically achieved by considering a super-evolution equation 

= , (2.28) 

where V is the supertime derivative with the defining property V 2 = —id/dt and Q is 
the supercharge of the theory satisfying the relation VQ + QV = 0. Then, the super- 
evolution equation provides a consistent square root of the Schrodinger equation, because 
it implies 

-i— = VH = VQm = -QVm = -Q 2 ^ = -Hm , (2.29) 
ot 

where H = Q 2 is the Hamiltonian of the system. Of course, the order of the equations 
that govern geometric evolutions is different, and there are also non-linearities that make 
the problem at hand to be more complicated mathematically. 

The analogy will be made more precise later by introducing anti-commuting parame- 
ters into the theory of geometric flows which promote time derivatives to super-evolution 
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operators. This procedure, which works well in two dimensions, has far reaching conse- 
quences, as it will allow to cast the Calabi flow into zero curvature form and deduce its 
properties from those of the Ricci flow. It is for this reason that both flows are discussed 
in parallel while uncovering their algebraic structures in two dimensions. In the process, 
we will also be able to construct a hierarchy of higher order geometric evolution equa- 
tions, together with a hierarchy of new infinite dimensional Lie algebras, which will be 
used for their zero curvature formulation. The general method works in the same way 
for higher dimensional Kahler manifolds, in that there is a squaring pattern for the time 
derivative operator that enables to connect consecutive flows, as above, but there is no 
zero curvature formulation for them above two dimensions. Thus, most of our work will 
be subsequently confined to two-dimensional spaces and postpone until the very end the 
derivation of some more general results that are valid in all dimensions. 

2.5 Landing in two dimensions 

Next, specializing to two dimensions and using a system of conformally flat (Kahler) 
coordinates, 

ds 2 t = 2e^ z '^dzdz , (2.30) 

we find that the only non-vanishing component of the Ricci curvature tensor is R zZ = 
—dd&. Then, the two different geometric deformations take the following neat form for 
the conformal factor of the metric, &(z,z;t), respectively, 

Ricci flow : d t $ = A$ , (2.31) 

Calabi flow : <9 t $ = -AA$ . (2.32) 

The symbol A stands for the Laplace-Beltrami operator, with the analyst's sign, 

A = e'^dd , (2.33) 

which depends on <3> and introduces non-linearities into the system. These equations can 
be be defined on any Riemann surface of genus g, but for all practical purposes we can 
concentrate to spaces with spherical topology, S* 2 . 

The dissipative character of these equations can be inferred from their linearization 
in the weak field approximation, where the curvature perturbations are assumed to be 
small. The reference metric around which such analysis should be carried, in order 
to make good physical and mathematical sense, depends on the flow. For the Ricci 
flow, it is appropriate to consider small fluctuations around the flat two-dimensional 
plane characterized by small and slowly varying values of $. This region corresponds 
to R^y ~ 0, that is to the asymptotically free limit of a non-linear sigma model with 
target metric g^ v . In the language of renormalization group flows, it corresponds to the 
ultra-violet region when the curvature of the space is positive and to the infra-red region 
when the space is negatively curved. Then, since $ ~ 0, the Ricci flow equation can 
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be approximated by the linear heat equation 9 t $ = <9<9$ in two-dimensional flat space, 
whose dissipative properties are well known. 

On the other hand, for the Calabi flow, it is more appropriate to consider the round 
sphere S 2 as reference frame for introducing small curvature perturbations around it. The 
choice is essentially uniquely and dictated by the limiting extremal metric of the two- 
dimensional Calabi flow and, therefore, one is led to examine fluctuations of the metric 
with g — (1 + /)<7o> where go corresponds to the line element ds 2 = 2Rf ] dzdz / '(1 + zz) 2 
of a sphere of radius Rq and / is taken small. Let us consider, without great loss of 
generality, axially symmetric deformations of the round sphere, which are parametrized 
by Legendre polynomials Pj(£) in the system of spherical coordinate (0,0), 

ds 2 t = R 2 [l + ei(t) PticosO)}^ 2 + sin 2 6d(f) 2 ) , (2.34) 

with corresponding small parameters e\{t). Such configurations qualify as approximate 
solutions to the Calabi flow equation, as they preserve the volume to linear order in ei, 

K = / dV(g e ) = 4nR 2 + 2nR 2 e l (t)(-l) 1 J* d£ = 4nR 2 (2.35) 

for all I > 1. In fact, the cases I = and I = 1 will be omitted in the calculations, as 
they both represent round spheres. Thus, we assume that I > 2. 

Then, to linear order in the perturbation parameters, the Calabi flow yields the 
following evolution 

e,(f) = c,(0)exp (-^*(* 2 -m + 2)) • (2.36) 

The result, which approximates well the asymptotic behavior of the full non-linear evolu- 
tion equation, when t — > +oo, shows that all perturbations are damped exponentially fast 
and the configuration settles down to that of a round sphere. More general deformations 
corresponding to spherical harmonics Y" z m (0, 0) with m^O can also be considered, but 
the results are essentially unaltered. All linearized solutions admit interesting physical 
interpretation in the context of multi-pole gravitational radiation from bounded sources, 
where they were first derived in the literature, as will be seen later in section 6. 

For comparison, a similar analysis is performed for the normalized Ricci flow, which 
is volume preserving and in two dimensions assumes the form 

d t $ = A$ + ^ (2-37) 

Kq 

when the space has spherical topology and its volume is taken equal to AtiRq. In this 
case, one considers linearized perturbations about the round sphere of radius R , which 
constitutes the equilibrium state of the normalized flow. Using the same ansatz (2.34) 
as before, it turns out that the corresponding parameters evolve, to linear order, as 

e,(t) = c(0)exp (-2^j(* " W + 2 )) > ( 2 - 38 ) 
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and therefore, as t — > +00, the perturbations are also damped exponentially fast 1 . Com- 
parison with the Calabi flow shows the same qualitative behavior, although in this case 
perturbations of a given angular momentum I diminish at a slower pace; the dependence 
of 6i(t) on /, being quadratic or quartic form, is a reflection of the corresponding evolution 
equations, being second and fourth order, respectively. 

It is natural to expect that the dissipative properties of the linearized flows will be also 
shared by the full non-linear parabolic equations of second and fourth order, respectively, 
as they do. The difficulty with the non-linearities is that they prevent the construction 
of explicit solutions, apart from a few isolated simple examples. In most cases, it is quite 
difficult to introduce ansatz that effectively reduce the non-linear equations to a simpler 
system of evolution equations for a finite number of moduli. The Ricci flow appears to 
be more tractable by such methods, since there is already a good number of simple (yet 
interesting) exact solutions. The situation is more primitive for the Calabi flow as only 
a handful of exact solutions is known. 

3 Algebraic framework 

Our primary aim in the following is to device a suitable algebraic method for establishing 
the integrability of the two different geometric deformations in two dimensions. It proves 
useful to sketch first the algebraic techniques that have already been used successfully 
for the formal integration of the Ricci flow and then extend their domain of applicability 
to the Calabi flow and further beyond. 

The notion of continual Lie algebras provides the main framework for the systematic 
formulation of the geometric flows as Toda field equations in two dimensions. The key 
idea is to treat the deformation variable t as a continuous "Dynkin index" of an appro- 
priately chosen infinite dimensional algebra and rewrite the Ricci and Calabi flows as 
zero curvature conditions in the two-dimensional system of conformally flat coordinates. 
Here, we present the basic definitions of continual Lie algebras and the Toda systems 
associated to them. We also review the general algorithm for obtaining their general 
solutions in terms of two dimensional free fields, using group theoretical techniques for 
the Backlund transformations. The specific choice of continual algebras that enable to 
cast different geometric flows into zero curvature form will be determined in the following 
sections together with the formal power series solutions of the equations. 

1 Of course, the time of the normalized flow depends logarithmically on the time of the unnormalized 
flow, as — i?olog(— t), according to the general relation between the two; thus, a spherical configuration 
that starts from very large size in the ultra-violet region t — > — 00 and diminishes linearly to zero size 
at t — by the Ricci flow, will correspond to a volume preserving deformation of the sphere under the 
normalized Ricci flow that now runs from —00 to +00 and tends to become round exponentially fast. 
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3.1 Continual Lie algebras 



The basic theory of continual Lie algebras is formulated by introducing a system of 
Cartan-Weyl generators {H(t), X ± (t)} that depend on a continuous variable t and satisfy 
the system of commutation relations, [15, 16], 

[X + (t),X~(t>)] = S(t,t')H(t') , [H(t),H(t')} = , 
[H(t),X ± (1?)]=±K(tJ)X ± (1>) . (3.1) 

K(t, t') is called the Cartan kernel of the algebra and it generalizes the Cartan matrix 
of simple Lie algebras to the continuous case, whereas the typical choice for S is S(t, t') = 
S(t—t'), which corresponds to 5ij in the case of discrete roots. This algebraic structure can 
be alternatively defined by smearing the Cartan-Weyl generators with arbitrary functions 
ip(t) of compact support, as in the theory of distributions. Thus, setting 

A(<p) = J dt(p(t)A(t) (3.2) 

the basic commutation relations assume the form 

[X + ^),X-^)} =H(SM)) , [H{ip),H{rj})\ =0 , 
[H{ip),X ± {^)\ = ±X ± (K(tp, ^)) (3.3) 

In general, the operators K and S are thought to be bilinear maps on the vector space 
of functions ip(t) satisfying the relations 

K(<p, Kty, x)) = Kty, K(<P, X)) , S(<p, Kty, x)) = S(K^, tp), X ) (3.4) 

for consistency with the Jacobi identities. Here, we only need to consider the case 
K(ip,ip) = (Kip)ip and S((p,ip) = S(ipip) prescribed by linear operators K and S, which 
always ensure the validity of the Jacobi identities. The corresponding continual Lie al- 
gebra will be denoted by Q(K,S), but it can be easily shown that Q(K,S) ~ Q(K,1) 
with K = KS; for it suffices to define a new basis of Cartan generators H(tp) = H(Stp). 
Therefore, for all practical purposes, the continual Lie algebras can be brought into stan- 
dard form, which is abbreviated to Q(K), setting S = 1 and dropping the tilde from K. 
From now on we will use this notation unless stated otherwise. The Cartan kernel K can 
be symmetrizable, but for the purposes of the present work we will also allow for algebras 
with anti-symmetric Cartan kernels K(t, t') = —K(t', t). In fact, the geometric flows turn 
out to be naturally related to continual Lie algebras with anti-symmetric kernels, as will 
be shown later in detail. 

The continual Lie algebras Q{K) are infinite dimensional Z-graded algebras with 
direct decomposition into subspaces 

Q(K) = (BnezGn ■ (3.5) 
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The Cartan-Weyl generators H and X span the components Go and G±i, respectively, 
thus forming the local part of the algebra G-i ©Go ©G+i, whereas all other elements of it 
can be constructed recursively by taking successive commutators of the basic generators 
so that Q n = [G n -i,G+i] f° r all n > 1 and G n = [£/n+i,£/-i] for all n < —1. The growth 
of the algebra which is characterized by the dimension of the subspaces Gn relative to 
the dimension of Go can vary accordingly depending on the choice of K. The derivation 
of all commutation relations of the algebra G(K) is not an easy task as there is no 
known analogue of the Serre relations in the general case. As a result, the complete 
description of G(K) is still lacking for arbitrary K, and in particular for those cases that 
will be encountered in the algebraic description of geometric flows. Nevertheless, the 
zero curvature formulation of the flows, and their integration, will be solely based on the 
local part of the algebras without ever requiring knowledge of the full structure of G{K) 
in either case. 



3.2 Associated Toda systems 

The Toda system associated to any continual Lie algebra with Cartan operator K is 
defined to be 

dd<S>(z, z;t) = J dt'K(t', t)e* (2 ' 2_;t,) . (3.6) 

$(z, z; t) can be viewed as a one-parameter family of two dimensional fields that depend 
on a continuous variable t rather than on discrete indices i, [15, 16]. As such, they 
generalize the Toda field equations associated to simple Lie algebras with Cartan matrix 
to the case of infinite dimensional algebras with a continuous system of roots, at least 
formally. In analogy to ordinary Toda systems, this equation is cast into zero curvature 
form 

[8 + A+(z, z), 5 + A.(z, z)\ = (3.7) 

using appropriate gauge connections A± taking their values in the local part Go © G±i of 
the algebra G(K). 

More precisely, let us consider the general ansatz for the two-dimensional gauge con- 
nections 

A±(z,z) =H(u ± )+X ± (f ± ) , (3.8) 

which depend on functions u± and f± of z, z and t. Some of these functions are superflu- 
ous and they can be eliminated by employing the gauge invariance of the zero curvature 
condition under the transformations 

A + ^g-\d + A + )g , A_ - g~\d + AJ)g (3.9) 

with arbitrary group elements g = expif (ip). As it turns out, it is always possible to 
choose the gauge 

A+ = H(V)+ X + (l) , A_=X-(e*) (3.10) 
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and then use the commutation relations of the algebra in order to obtain the following 
system of equations from the zero curvature condition, 

5^ = e* , <9$ = K(V) . (3.11) 

Elimination of the function \1> leads to the continual Toda field equation 

<9<9$ = tf(e*) , (3.12) 

which is the smeared form of equation (3.6) above. 

The general solution is obtained in analogy to ordinary Toda systems (for a review, 
see, for instance, [17]), by introducing a normalized highest weight state \t > that depends 
on the continuous variable t, 

X + (t')\t>=0, <t\X'(t') = 0, H(t')\t>=5(t-t')\t> (3.13) 

with < t\t >= 1. Then, the one-parameter family of two dimensional free fields with the 
decomposition 

$ (z,z;t)=f(z;t) + f(z;t) (3.14) 

following from <9<9$ = 0, is used to express the Toda field configurations in the form, 
[15, 16], 

$(z,z;t) = $ (z,z;t) -K (log < t\M-\z)M_(z)\t >) , (3.15) 

where M± are given by the path-ordered exponentials of the Lie algebra elements in Q± 1 , 
respectively, 

M+(z) = Vexp(^J Z dz' j dt'e f{z '' t ' ) X + (t , )j , 

M_(z) = Vexp^J* dz' J dt'e^'^X-^y (3.16) 

This formula summarizes the Backlund transformations of continual Toda systems, thus 
providing a free field realization of their general solution. 

The details can be worked out by expanding M± in power series and compute the 
expectation value of their product in the highest weight state, as 



< t\M- l M^\t >= 1 + ^(-l) n / dz x --- j~ dz n dz!--- "~ dz n x 

71=1 J J J J 

II dU / I] dt' t expf( Zl ; t i )exp/(^; ]n )> ■ - ■>-> ■-> ^ (3.17) 
i=i J i=i 

with the ordering z > z± > • • • > z n _\ > z n and likewise for the z's. The structure of the 
algebra Q (K) is fully encoded into the elements 

D {HM,-,t n -AA^-A} =< t \ X + (t 1 )X + (t 2 ) ■ ■ -X+(t n )X-(0 • ■■X-(t' 2 )X-(t[)\t > (3.18) 
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that depend on the choice of Cartan operator K. They are determined using the recursive 
relations 



^{*i,*2,— >*2'*">*»} = £ 5{tn _ g 5{t _ g _ £ K{t ^ q ) D 

j=i \ k=i 



{ti,"-,t n _i;t^,---,t^,"-,t^} 
i 



(3.19) 

with D t 1 ' 1 = 5(t — ti)5(t — t[). Here, tj is used to denote that t'j has been omitted by 
contraction with t m while shifting X + 's to the right and X^'s to the left. 

Thus, straightforward extension of the group theoretical methods used for the inte- 
gration of Toda systems for simple Lie algebras allows to construct the general solution of 
continual Toda field equations and obtain a formal infinite power series expansion around 
free field configurations. In regions where the exponential self-interaction becomes neg- 
ligible one simply has <£> ~ <3> for the corresponding Toda fields. Such formal expansion 
may or may not converge in the strict mathematical sense, but there is no general way 
to decide on such issues at the moment. 

The simplest and mostly studied example of this kind is provided by the so called 
heavenly equation, which is the following non-linear equation for a field $(z, z; t), 

J^e* = -dd* . (3.20) 

This equation is not parabolic, unlike the equations associated to geometric evolutions, 
for it involves second derivatives in space and time. It arises in the context of four- 
dimensional hyper-Kahler manifolds that admit a non-triholomorphic isometry associ- 
ated to a rotational S 1 action, and describes the self-duality condition on the Riemann 
curvature tensor of the metric in adapted coordinates, [18]. In the present context, it is 
identified with the Toda field equation for the large N limit of the SU(N) algebra, where 
the Cartan matrix = 25ij — 5i + ij — 5ij + i is replaced by the kernel K(t, t') = —5"(t—t') 
using a continuous parameter t to label the vertices of the Dynkin diagram, [19]. Thus, 
the corresponding Cartan operator is K = —d 2 /dt 2 . Integration of the heavenly equa- 
tion by group theoretical methods has been carried out, as for Toda theories, and various 
geometrically interesting solutions have been studied explicitly; the examples include the 
Eguchi-Hanson, Taub-NUT and Atiyah-Hitchin metrics, which all admit at least one 
non-triholomorphic S 1 action. In that context, the validity of the formal power series 
expansion of the general solution in terms of free fields has been tested in several cases, 
and the interested reader is referred to the literature, [20], for details. 

3.3 Generalization to super continual Lie algebras 

Next, we introduce the more general class of super continual Lie algebras which are defined 
by the basic system of commutation relations 

[X + {F),X-(g)} = E{TQ) , [E{T\E{Q)\ = , 

[H(F),X ± (g)] = ±X ± ((KF)Q) (3.21) 
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for the generators H and X^. The difference with the ordinary case is that the smearing 
functions T are taken to be super-functions on i? 1 ' 1 superspace with coordinate T = (t, 6) 
having 9 2 = 0. Likewise, the generators H(T) and X ± (T) are not bosonic but they also 
admit an expansion in it! 1 ! 1 superspace. However, the smeared form of the generators 
H(F) and X ± {T) is considered to be bosonic so that they possess commutation rela- 
tions rather than anti-commutation relations. Thus, in effect, supercontinual Lie algebra 
provide a special class of infinite dimensional Lie algebras that can be brought into a 
standard Cartan form with the aid of a Grassmann variable 9. 

The corresponding Cartan operator K can be an ordinary linear operator which is 
defined on the space of functions of a continuous variable t and acts separately on each 
component of the smearing super-functions T . More generally, it can be a super-operator 
acting linearly on the space of super-functions JF(T). In either case, the Cartan operator 
has to satisfy the relations 

[KT)Q = Q{KT) , (KT){KQ) = [KQ)[KJ=) (3.22) 

for all JF and Q, as follow from the Jacobi identity for the two triplets (H, X + ,X~) 
and (H, H, X ± ). These put no special restrictions on the choices that are available, 
and, therefore, the class of super-continual Lie algebras is well defined for all K. A 
particularly interesting example, which is relevant for the algebraic description of the 
Calabi flow, corresponds to the odd choice K = d/d6 — 88/ dt, as will be seen later. 
Other choices are also possible, but their physical relevance is not known in general. 

The smearing super-functions T can be expanded in components as 

F(T)=Mt) + 0Vi(t) , (3.23) 

using the superspace coordinates T = (t,6). Likewise, the generators of the algebra can 
be viewed as being parametrized by a continuous super-time parameter T with respect 
to which they also admit expansion into components, 

A{T) = A (t) + OAiif) , (3.24) 

with A being either H or X ± . By definition, the two components (fo(t) and y?i(t) are 
even functions taking values in ordinary numbers, and similarly Ao(t) and Ai(t) are also 
even. The smearing is defined by complete integration in i? 1 ' 1 superspace with measure 
dT = dtd9, as follows, 

A{T) = J A{T)T{T)dT = A ( Vl ) + A 1 ( l p Q ) , (3.25) 

using the known identities / d6 = and / 6d6 = 1. Thus, all such A(jF) are even, 
as required. Then, the structure of the commutation relations can be worked out in 
components for the bosonic generators Hi((p) and Xf(tp) with % — 0, 1, depending on the 
form of K. 

Toda systems associated to such algebraic structures can be defined as usual by the 
zero curvature form [d + A + (z, z), 8 + A_(z, z)\ = 0, choosing 

A + = H{Q) + X + (l) , A_=X~(e^) (3.26) 
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for the gauge connections. Then, if the super-functions JF and Q are related to each other 
as dT = K(Q), the corresponding Toda field equation will be 

ddT = K{e r ) . (3.27) 

Expanding JF in components and comparing the even and odd parts of the equation 
yields an equivalent description of the associated Toda system. The general solution can 
also be obtained in this case by straightforward generalization of the group theoretical 
methods developed for ordinary Toda systems. Further details will be given later for the 
particular example of the Calabi flow. 

Other possibilities may arise, in principle, by considering JF as super-functions on 
R 1][N superspace, according to our general discussion. The structure of the associated 
supercontinual Lie algebras can also be worked out in component form, but become 
more complicated for higher vales of N. 



4 Ricci flow 

According to the general discussion in section 2, the Ricci flow equation in two dimensions 
assumes the following simple form in a system of conformally flat coordinates (z, z), 

dd${z, z- 1) = $e* ( *- 2;t) , (4.1) 

where $ is the conformal factor of the metric undergoing continuous deformations. As 
such, it admits a natural algebraic interpretation in the framework of Toda field equations 
for the choice of continual Lie algebra with Cartan operator K = d/dt. In particular, 
following [21, 22], we consider the continual Lie algebra whose local part satisfies the 
commutation relations, in smeared form, 

[X + (y?), X-ty)] = H(^) , Hty)] = , 

[H( V ),X±m = ±X±(^) ■ (4.2) 

Then, the two-dimensional Ricci flow, which is viewed as Toda system for this algebra, 
[21, 22], admits the zero curvature formulation 

[d + A + (z, z), 5 + A.(z, z)\ = , (4.3) 

where A± take values in the local part as 

A + = H^)+X + (l) , A_ = X"(e*). (4.4) 

Also, following the general discussion in section 3, we find that (9$ = d t ^ and = exp$ 
and, hence, $ satisfies the Ricci flow equation (4.1) after eliminating ^. 

The main ingredient in this particular formulation of the Ricci flow is the use of a 
continual Lie algebra to incorporate the deformation variable t into its root system. Thus, 
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the equation is integrable in two dimensions using flat space coordinates (z,z), without 
reaching a contradiction of terms with the dissipative behavior in time. In this context, 
the parabolic nature of the flow, which is first order in t, is held responsible for the novel 
structure of the underlying Lie algebra with anti-symmetric Cartan kernel 

K(t,t') = -p(t-t>). (4.5) 

Such algebras have been studied very little in the mathematics literature, but it has 
already been noted that Q(d/dt) exhibits exponentially fast growth beyond its local 
part. As a result, its complete structure is still lacking and there is no framework in 
which to realize transformation rules for all the generators. Also, the precise meaning of 
integrability for systems associated to such algebras should be understood better in the 
future, and in connection with the explicit construction of their conservation laws in two- 
dimensional space. The fact is that they become relevant for the algebraic description 
of certain dynamical problems with the Ricci flow being the first example of this kind. 
The Calabi flow also falls in this category for appropriate choice of Cartan operator in 
the class of supercontinual Lie algebras, as will be seen later in detail. 

Next, we apply the general formalism to describe the solutions of the two-dimensional 
Ricci flow equation in terms of a one-parameter family of free fields $o(z, z; t) = f(z; t) + 
f(z;t). The general expression (3.15) specializes to 

$(z,z;t) = %(z,z;t) - d t (log < t\M?(z)M_{z)\t >) , (4.6) 

where M± are the path ordered exponentials (3.16) that should be evaluated for the 
particular algebra Q{d/dt). This provides the solution in closed form by group theoretical 
methods, but it is very implicit. In practice, one should evaluate recursively the elements 
D^ tl,t2 ' •* n, *i' i 2' »*«} an j substitute them into the general power series expansion (3.17). 

Explicit calculation yields in this case the following results, up to n = 3, [21], 
D? 1 ''^ =S(t-t 1 )S(t 1 -t[) , 

D {tutrAA) = 5{t _ ti)5[ti _ t r )5{t2 _ q {25{t _ h) _ dt5[t _ h)) 
D {t M ,t'M} = s{f _ fi)§{ti _ t r )6{t _ h)5[t2 _ Q §{h _ q x 

x (28(t-t 2 )-d t 5(t-t 2 )) + 
+ 8(t - hWh - t[)5(t 2 - t' 3 )8(t 3 - Q x 

x (S(t - 1 3 ) - d t S(t - t 3 )) (2S(t - 1 2 ) - d t S(t - t 2 )) + 
+ 5(t - *i)<S(ti - t[)5(t 2 - t' 2 )5(t 3 - 1' 3 ) x 

x (25(t - t 2 ) - d t S(t - 1 2 )) (S(t - t 3 ) - d t S(t - t 3 ) - d t2 S(t 2 - t 3 )) (4.7) 

and so on for higher n. The expressions become considerably longer beyond n = 3 and 
they are not included here. 
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In special cases, where the corresponding geometric deformations are axially sym- 
metric, the solutions can be taken to depend only on the combination of coordinates 
z + z := Y. Then, for free field configurations of the form $ (Y;t) = cY + d(t), the 
systematic expansion of the path ordered exponential leads to the power series solution 

$ = $o + + (2!^* (^*°) + j^dt (3e*°(d t e*°f + e 2 *°9 2 e*°) 

+ W^f 9t ( e ^° d t^° + 17e 2 *°(d t e*°)(c> 2 e*°) + 18e*°(c> t e*°) 3 ) 

+ (5!^j2' 9 * ( e4 *° 9 * 4e *° + 36e 3 *°(9 t e*°)(9 t 3 e $0 ) + 35e 3 *°(d 2 e*°) 2 

+ 324e 2 * (9 t e*°) 2 (9 2 e* ) + 180e*°(9 t e* ) 4 ) + • • • (4.8) 

by evaluating the contribution of all terms in the expansion up to order n = 5 included. 
It should be understood as a power series expansion of exp$ around the free field con- 
figuration exp$ that becomes very small when Y — > ±oo for c < or > 0, respectively. 
It can also be verified directly that this expansion satisfies the Ricci flow equation order 
by order in powers of exp$ . 

The validity of the formal power series expansion and its convergence have be tested 
in various examples of axially symmetric geometric deformation that are known in closed 
form using more direct (mini-superspace) methods. A notable example is the sausage 
model that corresponds to the choice of conformal factor, [23], 

e* (y;t) = — r-^T — - with < X < 2tt , - oo < F < +oo , (4.9) 

a(t) + 6(t)cosh2F ~ y J 

while t is running from — oo to and 

a{t) = - 7 coth(2 7( ) , Kl) = . (4.10) 

Note that a > b > since t assumes only negative values; t — > — oo corresponds to the 
asymptotically free ultra-violet limit of the 0(3) sigma model in two dimensions. This 
solution describes a family of axially symmetric deformations of the round 2-sphere for 
all values of the parameter 7 > 0. It fits the expansion (4.8) in power of exp(— 2Y) for 
large Y, [21], provided that the free field is chosen to be &o(X', t) = cY + d(t) with 

c = -2 , d(t) = log ^-^sinh(2 7 t)^ . (4.11) 

Thus, there is mounting confidence that the formal manipulations used for the inte- 
gration of the Ricci flow indeed make good sense. Further details can be found in the 
published works, [21, 22]. The only drawback of the present formalism is the inability to 
prescribe directly the geometric shapes of space from different choices of free fields. 
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5 Calabi flow 



Next, we turn to the two-dimensional Calabi flow, which will be shown to admit zero cur- 
vature formulation for appropriately chosen gauge connections. The method of investiga- 
tion involves the introduction of anti-commuting variables following a simple observation 
that formally relates it to the Ricci flow. Then, we give full justification for the choice of 
the infinite dimensional algebra that enables to cast the Calabi flow in integrable form, as 
Toda system, and extend the validity of the group theoretical methods introduced earlier 
in order to describe (at least formally) the general solution in terms of free fields, as for 
the Ricci flow. The results described below complement earlier work on the subject by 
providing full justification of the seemingly ad hoc choices made in reference [24] and go 
far beyond it in some directions. 



5.1 Calabi flow as super-evolution 

According to the general discussion found in section 2, the two-dimensional Calabi flow 
assumes explicitly the form 

d t e* = -dd (e-*dd<S>) (5.1) 

in a system of conformally flat coordinates with conformal factor <&(z, z; t) varying with 
time t. Specializing an earlier observation, it is instructive to compared it to the two- 
dimensional Ricci flow equation (4.1), which is equivalently stated as 

^e* = dd (e-*dd<S>) (5.2) 

after taking its derivative with respect to t. Indeed, the Calabi flow appears to be formally 
related to the Ricci flow by taking the square root of its evolution operator d/dt, up to 
a sign. 

This simple correspondence motivates us to introduce an anti-commuting variable 9, 
with 9 2 = 0, as the supersymmetric partner of the deformation variable t and write down 
the following abstract equation 

ddT = V T expF (5.3) 

where 

T{z, z; t, 9) = $(z, z; t) + 9V(z, z; t) (5.4) 
is a mixed superfield with bosonic components $, \& and 

v - = m ~ e wt < 5 - 5 > 

is the associated super-derivative in i? 1 ' 1 superspace with coordinates T = (t,9). It 
satisfies the relation T>\ = —d/dt, whereas the complementary super-derivative operator 
d/d9 + 9d/dt squares to d/dt, since it relates t to —t. The equation above can be 
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viewed as some short of super-evolution equation of Ricci type with respect to the super- 
time T. Its precise content can be extracted by working in components. Expanding 
expjF = (1 + #^)exp$ and comparing the even and odd parts, it readily follows that 

e*^ = dd$ , «9 ie * = -dd^ . (5.6) 

Subsequently, by eliminating it follows that <£> satisfies the Calabi flow equation (5.1) 
above. Thus, by taking the square root of the time derivative operator allows to connect 
the two distinct classes of geometric deformations of second and fourth, respectively, via 
super-evolution, as was anticipated before on general grounds. 

This simple result has far reaching consequences, as it allows to prove integrability 
of the two-dimensional Calabi flow by viewing it as Toda field equation for a particular 
supercontinual algebra with odd Cartan operator K = T>t, in analogy to the Ricci flow 
that corresponds to the choice K = d/dt. More precisely, we are led to introduce an 
algebra whose local part satisfies the commutation relations 

[X + (F),X-(g)] = H{TQ) , [H(f),H(g)] = , 

[H^X^G)] = ±X ± ((V T F)Q) . (5.7) 

It defines a supercontinual Lie algebra in the nomenclature introduced earlier. 

Then, following the general discussion in section 3, we write down the zero curvature 
condition 

[d + A+(z, z), 8 + A-(z, z)] = 0, (5.8) 
where the connections A± take values in the local part of the algebra, as 

A + = H{G) + X+(l) , A_ = X-(expJ^) . (5.9) 

The superfield T satisfies the flow equation (5.3), provided that Q is constrained to 

dT = V T Q (5.10) 

by the zero curvature condition. Thus, the Calabi flow is indeed an integrable two- 
dimensional system of Toda type. 

Next, the structure of the algebra (5.7) is analyzed in more detail, using components 
in i? 1 ' 1 superspace, and then used further to integrate the equation by group theoretical 
methods, as in ordinary Toda systems. 

5.2 Working in components 

The algebra introduced for the Calabi flow is not a supersymmetric extension of what 
was used before for the Ricci flow, but it is bosonic. Our prescription did not involve any 
substitution of the space derivatives d and 8 by the supersymmetric covariant derivatives 
T>±, respectively, as it is customary in supersymmetric systems, but only used the odd 
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generalization of the time evolution operator in i? 1 ' 1 superspace, in order to converge 
a second order differential equation in (z, z) into fourth order. It can be alternatively 
viewed as changing the kernels of the algebra for the Ricci flow, 

S(t,t') = 5(t-t') , K{t,t') = -d t 5(t-t') , (5.11) 

into the two super-kernels, respectively, 

S{T, T') = 5{T - T') = (6- 6')5(t - t') , (5.12) 

K(T, T') = -V T 5(T - T') = -5{t - t') - 69'd t 5(t - t') , (5.13) 

thus giving rise to the bosonic super continual Lie algebra with K = T>t- Of course, the 
passage from smeared to unsmeared form of the algebra has to be done with care, since the 
integration over the Grassmann coordinates 6 and 6' will flip signs in some places. This 
generalization also changes the symmetry of the corresponding kernels, since S(T, T') is 
now anti-symmetric and K(T, T') is symmetric under the interchange of their arguments, 
T <-> T. 

It is useful in many respects to work out the form of the algebra in components by 
expanding the generators, 

X ± (T) = X±(t) + exf{t) , H(T) = H (t) + 0ifi(f) . (5.14) 

All components Xf(t) and H^t) with % — 0, 1 are bosonic, and one arrives to the equiv- 
alent description of the local part of the algebra with commutation relations 

[H (t), X±(t')\ = ±5(t - t')X±(t>) , [X±(t), Xf(t')\ = ±5(t - t')H (t>) , 

[H (t), Xf(t')} = ±6(t - t')X±(t') , [X+(t), Xf(f )] = S(t - t'^i?) , 

[H^t), X±(t')} = ±d t 5(t - t')X±(t') , (5.15) 

whereas the remaining ones are trivial, i.e., 

[H^t), X±(t>)] = , [X + (t), x -(O] = , 

[Hi(t), H J (t')] = 0. (5.16) 

Smearing with ordinary functions of the continuous variable t, the algebra takes the 
form 

[H (<p), XtW] = ±X?(<rt) , [X+(<p), X-{i>)\ = H^) , 



and 



[H^ip), X ± (^)] = , [X + (y?), X^m = , 

[Hi(<p), Hjty)] = . (5.18) 
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One could have started directly from such an algebra, without prior justification, as 
in reference [24], and use it to establish the integrability of the Calabi flow 2 . However, 
employing the notion of super continual Lie algebras proves useful for putting its structure 
in neat form and making direct connection with the special class of integrable systems 
of Toda type. 

The zero curvature formulation of the Calabi flow based on the algebra (5.7) can also 
be analyzed in components. The gauge connections (5.9) take the form 

A + = H (f) + H^g) + X+(l) , A_ = X (Ve*) + X~(e*) , (5.19) 

using the expansion of the mixed superfields T = $ + 0\& and Q — g + Of. Also, the 
relation (5.10) among these fields reads in components as 

d$ = f , d^ = -d t g . (5.20) 

As a result, the zero curvature condition implies two more equations that follow by 
comparing coefficients of the terms H and Hi, respectively, 

Bf = ^e* , Bg = e* , (5.21) 

and lead to the system (5.6) that accounts for the Calabi flow, after eliminating all other 
functions but <E>. 



5.3 Integration of Calabi flow 

The general solution can be obtained in analogy to the Ricci flow equation, as in ordinary 
Toda systems, by introducing a normalized highest weight state |T >, assuming that it 
exists, which depends on the continuous coordinate T of i? 1 ' 1 superspace with < T\T >= 
1. Equivalently, in terms of components (t,9), we may assume that the highest weight 
state is expanded as 

|T >= \t > +9\t >i (5.22) 

with the orthonormalization conditions that the bosonic state \t > has unit norm and 
it is perpendicular to \t >i. It is more convenient to present its defining properties in 
smeared form 3 , 

x+(r)\g>=o, <g\x-(F) = o, h{f)\q >= \j=q > , (5.23) 

2 There is a slight difference of signs in the commutation relations of the algebra and those appearing 
in reference [24]. The two choices are simply related to each other by changing Xq — ► and 
Xf — ► =fiX^ and they do not affect the proof of integrability. The signs that arise here guarantee that 
the algebra is compatible with the Hcrmiticity condition h\ = Hi and X^ = Xf for alH = 0, 1, unlike 
the other case. These Hcrmiticity conditions turn out to be essential for the actual integration of the 
equation by group theoretical methods using highest weight representations, as will be seen later. 

3 Note that we systematically avoid to write the commutation relations of the supercontinual Lie 
algebra and the definition of its highest weight state starting from the unsmeared form of the super- 
generators A(T) labeled by T in i? 1 ' 1 . In this way, various sign discrepancies associated to the ambiguity 
of ordering products of super-functions defined at different points of super-space are avoided, as they 
can easily lead to errors. 
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where A(T) = A Q {ifi) + Ai{<p Q ) for the generators of the algebra and 

\T >= J\T> F(T)dT = J (y>i(f)|f > +<p (t)\t >i) dt = \ Vl > + |y? >i (5.24) 

for the highest weight state, using the smearing super-function JF(T) = <po(t) + 9ipi(t). 
The components of equation (5.23) satisfy the smeared relations 

H ((p)\il> > = , H (<p)\ip >i= \(pip > , 

H^Hj >i= >i for i = 0, 1 , (5.25) 

and 

X+( V )\^> t =0 = X+( V )\^> i for 2 = 0,1 (5.26) 

together with their Hermitian conjugate equations, using ordinary functions ip(t) and 
ip(t). The same relations can also be written in unsmeared bosonic form, as 

if (OI*>o=0, H (t')\t> 1 =5(t-t')\t> , 

H l (t')\t> i = 5(t-t')\t>i fori = 0,1, (5.27) 

and 

X+{t')\t >i= = X+(t')\t >i for % = 0, 1 . (5.28) 

It follows immediately from above that H (t')\T >= 95 (t — t')\t > and Hi(t')\T >= 
5(t — t')\T >, which will be useful later for the calculations. These relations also imply, in 
particular, the following condition for the unsmeared action of the Cartan generator on 
the highest weight state H(T')\T >= {9 + 9')5{t-t')\T >, which is neither 5{T-T')\T > 
nor |T > 5(T — T') as compared to naive expectation. 

The zero curvature formulation of the Calabi flow, which is better stated as super- 
evolution equation (5.3), allows us to express its general solution in superfield form 

T(z, z; t, 9) = F (z, z; t) - V T (log < T\M-\z)M_{z)\T >) , (5.29) 

where JF = <3> + 9^ Q is a one-parameter family of mixed superfields whose components 
are two-dimensional bosonic free fields that depend on the continuous variable t, 

$ (z,z;t) = f(z;t)+f(z;t), V (z,z;t)=4>(z;t)+4>(z;t). (5.30) 

Also, setting Tq{z\ T) = f(z; t) + 9ip(z\ t) and T${z\ t) = f(z; t) + 9tjj(z; t), the operators 
M± are given by the following path-ordered exponentials 

M + (z) = Vexp^J* dz> J e^ z '' T ')X + (T')dT>y 

M_(z) = Vexp^J 2 ' dz' j 'e^ s '^X-(T')dT'y (5.31) 
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Performing the integration over the anti-commuting variable 9' that enters in the 
exponent, one arrives at the equivalent expressions in terms of components, respectively, 

M+(z) = Vexp^ dz' J dt'e fiz '' tr) [X+(t') +ij(z'-t')X+(t')]y 

M_(z) = Vexp (J* dz' J dt'e^'^iX^it') + i')* ~(O]) . (5.32) 

Then, formula (5.29) summarizes the Backlund transformation of the Calabi flow, as a 
continual Toda system, and provides realization of its general solution for the arbitrary 
free field configuration ($ , ^o)- The solution for $(z, z; t) is obtained by comparing the 
even terms of the mixed superfield equation (5.29), whereas \J/(z, z; t) is obtained from the 
odd terms; the resulting expressions, which are only implicitly defined here, should satisfy 
the equation \I/ = exp(— $)<9<9$, by construction. The (formal) power series expansion of 
the solution follows by expanding the path-ordered exponentials, as usual. 

In the following we will consider first, for simplicity, the class of solutions with 
^o(z, z; t) = 0, as they also resemble the solutions of the Ricci flow equation in some 
respect. The possibility to include free field configurations with \& 7^ will be consid- 
ered separately, at the end of this section, as they lead to more complicated expressions 
for the free field expansion of the general solution. 

(i) Solutions with ^ = 0: In this case, the path-ordered exponentials (5.32) 
simplify considerably, as they involve only the generators Xf{t') in their exponent. Then, 
the detailed structure of the solution can be worked out in power series, which assumes 
the special form 

°° fZ fZn-l fZ fZn-l 

< T\M- X M_\T >= 1 + ^(-l) n / dz x ---\ dz n dz x ---\ dz n x 

n=l J J J J 

/n „ n 

]J dti / 1] dt^pf( Zi ; tOexp/fo fyDp'"'^'"'™ (5.33) 
i=i J i=i 

with the usual path-ordering prescription z > Z\ > ■ ■ ■ > z n _\ > z n , and likewise for the 
z's. The structure of the corresponding algebra (5.7) is fully encoded into the elements 

4 tl ' t2 '"'' tm;t ' 1 ' t ' 2 '"'' C} =< TlX+MX+fa) ■ ■■Xt{t n )X^{t' n ) ■ ■■Xr(t' 2 )Xr(t[)\T > 

(5.34) 

that only involve strings of the component operators X] t (t), for — 0. Of course, even 
in this case, the other components of the algebra generators, X^(t) as well as Hi(t), will 
also appear in the process of evaluating the corresponding matrix elements by pushing 
X + 's to the right and X~'s to the left, as they follow from the commutation relations. 

These elements can be computed recursively, as usual, starting from the defining 
relation < T\T >= 1. The result turns out to be 

Dp' A} = S(t - tJSih - t[) , 
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x (25(t-t 2 ) + 9d t 5(t-t 2 )) + 
+ S(t - t^Sih - t^Sih - t' 3 )S(t 3 - t' 2 ) x 

x (5{t - t 3 ) + 9d t S(t - t 3 )) {26{t - t 2 ) + 9d t S{t - t 2 )) + 
+ 5(t - ti)<J(*i - t[)5(t 2 - t' 2 )5(t 3 - t' 3 ) x 

x (25{t - t 2 ) + 9d t S(t - t 2 )) (S(t - t 3 ) + 9d t S(t - t 3 ) + 9d t2 5(t 2 - t 3 )) - 
- 5{t - t x )5{t - t 2 )5(t x - t[)8(t 2 - t' 2 )5(t 3 - t' 3 )d t 5(t - t 3 ) (5.35) 

and so on. Higher order terms with n > 3 are also calculable, but their expressions 
turn out to be quite lengthy with rapidly increasing complexity. For this reason they 
are not appended here, although contributions up to n = 4 and 5 will be included in the 
subsequent evaluation of the field $. 

The expressions above should also be compared to the matrix elements (4.7) corre- 
sponding to the Ricci flow, as they show many similarities for small values of n. Ob- 
serve that the time derivatives of delta functions appearing in the Ricci flow matrix 
elements now come multiplied with —9, which is consistent with the change of d/dt to 
T>t = d/d9 — 9d/dt as one passes to the algebra of the Calabi flow. However, there is 
also a new term appearing in the last line of the n = 3 matrix element, with no analogue 
in the Ricci flow, which contributes an additional term, without 9, involving the operator 
T>\ = —d/dt in the general power series expansion. The differences are becoming more 
visible in the structure of higher order terms, reflecting the more complicated nature of 
the commutation relations, which now involve a system of six rather than three basic 
bosonic generators of the algebra that all mix among themselves. 

It is relatively easy to perform the path-ordered integrations and work out the form of 
the power series expansion for axially symmetric deformations that depend only on the 
combination Y = z + z for which the free field has the simpler looking form &o(Y; t) = 
cY + d(t), as for the Ricci flow. For constant values of the parameter c, it turns the 
element < T\M+ l M_|T > can be computed as series expansion in powers of exp$ and 
the Toda superfield T assumes the following form 

T = % + j^Vre** + j^V T (e*°P T e*°) + -^V T (e™°V*e*°) 

+ 1 V T (e 3 *°^ e <i.„ + He 2<E, o(P Te *o)(p2 e <i. ^ 
{4. c ) 

+ - i- P T fe 4j, °P^ $ ° + 29e 3<E, °(p2 e *o)2\ + . . . (5>36) 
(5! c 5 ) 2 v 1 

by taking into account the contribution of the matrix elements (5.34) up to order n = 5. 
It is instructive to compare the terms with those appearing in the power series expansion 
(4.8) for axially symmetric deformations of the Ricci flow. To each order we see that the 
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terms transform from one into the other by simply changing d/dt into T>t- Polynomial 
terms of the form ((9 t exp$o) 2 > an d alike, are missing from the Calabi flow, since they 
will correspond to (Vxexp^ ) 2 and vanish for ordinary bosonic fields &o, due to 9 2 = 0. 
As for the remaining terms, there are some differences in their relative coefficients that 
appear to order n > 4 in the expansion and reflect the differences of the two algebras, as 
expected on general grounds. 

Substituting the result into the general expression (5.29) with \&o — an d comparing 
the coefficients of the even and odd terms, the expansion of $(Y; t) in terms of free fields 
reads 

$ = $o - J^^d^ + (4^452 (**°# e *° + H(^-) 2 ) + • • • , (5.37) 
whereas the expansion for ^(Y; t) reads, up to fifth order, 



* :: << = -7TT^ d ^° + (e*°c> t 2 e*° + 2{d t e^f) 



(1! 

( e 2*oa3 e *o + 62e*°(^e ,I>0 )( ( 9 t 2 e* ) + 87(9 t e*°) 3 ) + • • • . (5.38) 



5^2 



(5! c 

The validity of these expressions can also be verified directly, order by order, by substi- 
tuting into the system (5.6) that describes the Calabi flow equation. However, unlike the 
Ricci flow, there are no simple mini-superspace solutions that describe axially symmetric 
deformations of the round sphere by the Calabi flow in closed form and which resemble 
the sausage model. It remains an interesting open question to find such models, provided 
that there are consistent mini-superspace reductions of the Calabi flow, and compare the 
exact solutions with the formal power series expansion above. Similar remarks also apply 
to solutions with non-trivial free field component that are examined next. 

(ii) Solutions with ^ 7^ 0: In this case, the choice of fields should be compatible 
with the equations 

e*°* = dd% , dd^o = -— (5.39) 

when applied to the most general free field configuration ($ , \J/ )- Of course, both sides 
of these equations vanish, since <3> is a two-dimensional free field that approximates well 
$ only in a region of space where exp$ — 0, thus leaving undetermined. Such 
asymptotic behavior is commonly assumed in all theories of Toda type in order for the 
formal power series expansion of their solutions to be meaningful, if also convergent. A 
simple choice is to consider independent of the coordinates (z, z), as it will provide free 
field expansion for the configurations of Calabi flow around constant curvature spheres, 
with R = —^o, and in the vicinity of those points where $o — * — oo. The choice \l/o — 
that has been considered so far corresponds to such spheres, but with infinite radius, 
whereas in general one may assume that \l/o is non-zero, but it can also depend on time, 
^o(t). Even more generally, can be an arbitrary family of two-dimensional free fields 
that are independent of the choice $ - The analysis below is quite general, without 
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assuming any restrictions on the form of ^q, apart from the case of axially symmetric 
solutions that are only considered for illustrative purposes. Other interesting possibilities 
also arise by specializing to particular classes of free field configurations, as will be seen 
later. 

The calculations are considerably more involved now for evaluating the expectation 
value < T\MT l {z)M_{z)\T > term by term in the power series expansion. Already, the 
first non-trivial term arising at order n = 1 takes the form 

J dz[ J dz[ J dti J dt' 1 expf(z' 1 ;t 1 )expf(z' 1 ;t' 1 ) x 

x < TliX+ih) +^(z[;t 1 )X+(t 1 )][X^t> 1 ) +$(z[;t' 1 )X-(t' 1 )}\T >(5.40) 

as can be readily seen using the path-ordered exponentials (5.32) for M ± . For \17 ^ 0, 
there is an additional contribution coming from the cross terms in the expectation value, 
which were not present before. For higher values of n there are many more new terms 
that contribute to the final answer, making the process tedious. 

As before, it is also convenient here to restrict attention to axially symmetric configu- 
rations that depend only on the variable Y = z + z and assume that the two independent 
free fields are parametrized as follows, 

$ = c(z + z) + d(t) , = a(z + z) + b(t) , (5.41) 

where a and c are arbitrary constants whereas b(t) and d(t) depend generically on time. 
The path-ordered integrals can be evaluated easily in this case, but the resulting ex- 
pressions are more lengthy than the simpler case with \17 = 0. It is quite instructive, 
nevertheless, to present the contributions to order n = 1 and 2, 

< T\M?M.\T >= 1 - j^e*- (l + 2 -fj) + ^L_[e e *»a te *. 

- (»„ - | - 2) e* - B (*S " 4 (l + f) % + £ + e»°] , (542) 

whereas higher terms are of order exp(3$o) or more. 

Then, plugging into the general expression (5.29) and comparing the even and odd 
parts of the equation, it follows to this oder that 

(5.43) 

and 

*-*o-(T^^"-^((*»-T) e2t ")^ (5.44) 
It can be verified, as independent check, that these expansions are consistent, order by 
order in powers of exp$ ) with the dimensionally reduced equations <9y$ = ^/exp*^ and 
dty = — (9texp$ accounting for axially symmetric deformations along the Calabi flow. 
Higher order corrections can also be included, if required. 
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Finally, note in passing that there are special configurations in which both free fields 
$ and \l/o do not depend on time, t, and which can be used to provide power series 
solutions of the simpler equation AA$ = 0; equivalently they can be thought as describ- 
ing Backlund transformations of the non-linear equation <9<9$ = ^exp^ which looks 
like a generalized Liouville theory with space dependent coupling Our method is 
also capable to provide solutions of this equation in terms of an arbitrary free field <3>o 
for any given field \l/o- We will say more about this time independent class of configu- 
rations, and their formal power series solutions, while discussing the so called type III 
Robinson- Trautman metrics in four dimensional general relativity. 



6 Robinson- Trautman space-times 

In this section we discuss the embedding of the two-dimensional Calabi flow into four- 
dimensional vacuum Einstein equations and review the main characteristics of the space- 
time solutions that are known up to this day. This material is included to make the 
presentation self-contained, as it allows to put our results in the physical context of 
spherical gravitational waves. 



6.1 General structure of the metrics 

Recall that the most general gravitational vacuum solution in four space-time dimensions 
which admits a geodesic, shear-free, twist-free but diverging null congruence is given by 
the general class of Robinson- Trautman metrics, [25] 

ds 2 = 2r 2 e*dzdz - 2dtdr - Hdt 2 , (6.1) 

where H(z,z,r,t) has the special form 

H = rd t <$> - A$ - (6.2) 

expressed in terms of a single function $(z,z;t), and m(t) is a one-parameter family of 
parameters that in some cases represent the physical mass of the system. The affinc 
variable r varies along the rays of the repeated null eigenvector and t is a retarded time 
coordinate. Closed surfaces with constant r and t represent distorted two-dimensional 
spheres, and, therefore, t-dependent solutions are thought to describe spherical gravita- 
tional radiation; of course, these are not exactly spherical gravitational waves in vacuum, 
since spherical symmetry implies that the above metric is static. It is more appropriate 
to think of the Robinson- Trautman metrics as describing gravitational radiation outside 
some bounded region. 

There is a close analogy between Lienard-Wiechert fields of Maxwell equations and 
Robinson- Trautman metrics of Einstein equations in that they both admit a principal 
null vector field which is geodesic, shear and twist free with non-vanishing divergence, 
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thus giving a more intuitive meaning to the general ansatz that was made above, [26]. 
In the gravitational case, the field equations imply that all vacuum solutions in this class 
satisfy the following fourth-order differential equation for $(z, z; t), 

AA$ + 3md t $ + 2d t m(t) = , (6.3) 

where A = exp(— $)<9<9 as before. To understand its equivalence to the Calabi flow in two 
dimensions, (5.1), it suffices to relabel the null hypersurfaces of the Robinson- Trautman 
metrics and set m equal to constant whenever m(t) ^ 0. The physically interesting case 
corresponds to positive values of m, in which case t runs forward to infinity, but one 
may also consider the opposite case as well. Then, in the gauge with constant m, the 
Robinson- Trautman equation becomes 

AA$ + 3md t $ = (6.4) 

and it is identified with the Calabi flow by simple rescaling of the retarded time, so that 
3m = 1; as such it provides an exact embedding of the two-dimensional Calabi flow into 
four-dimensional gravity, as it was first noted in reference [27]. 

Finally, we note for completeness that there are also interesting generalizations of 
the present framework to Robinson- Trautman space-times with cosmological constant A, 
which are obtained by modifying the function (6.2) appearing in the metric to H\ = 
H — Ar 2 /3. Vacuum Einstein equations with cosmological constant reduce again to the 
same differential equation, as given in (6.3), which turns out to be independent of A. For 
an overview of the subject, see, for instance, [28]. 

6.2 Petrov classification of metrics 

Next, following [28], we review the Petrov classification of Robinson- Trautman vacuum 
solutions depending on the algebraic relations among the surviving components of the 
Weyl curvature tensor that are typically represented by the quantities 

* 2 = *3 = ^e-/ 2 a ( A$), 

^ 4 = (e-*d(rd t ® - A$)) (6.5) 



and their complex conjugate expressions, 
space-times according to the list 

Type N : 
Type III : 
Type D : 
Type II : 



There are four different algebraic types of 



# 2 = $ 3 = , 

^ 2 = 0, ^3^0, (6.6) 
3^4 = 2^ 3 2 , *2^0, 
3^4 ± 2^ 3 2 , ^2^0. 
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The mass parameter m, which is taken constant, vanishes for the first two types and as 
a result the Robinson- Trautman equation reduces to 



AA$ = , (6.7) 

with t appearing as spectator without apparent role in dynamics. The last two types are 
characterized by the condition m/0 and they correspond to genuine trajectories of the 
Calabi flow as given by equation (6.4). More precisely, we have the following: 

Type N solutions satisfy additional conditions, by their definition, which give rise to 
the special equation 

dd$ = K(t)e* , (6.8) 

using an arbitrary integration function of the variable t alone, K(t), that plays the role of 
coupling constant in the resulting two-dimensional Liouville equation. In this case, the t- 
dependence of space-times is not in general compatible with the Calabi flow, since m — 0, 
unless K(t) is constant. We only note here that flat space-time arises as special solution 
of this type, whereas non-flat type iV solutions typically have lines of singularities in the 
three-dimensional (r, z, z) space and they are also of limited interest in general relativity. 

Type III solutions satisfy the equation AA$ = 0, as before, with the condition 

<9(A$) ^ , (6.9) 

and its complex conjugate, by their definition. Therefore, space-times of this type corre- 
spond to solutions of the differential equation 

A$ = f(z;t) + f(z;t) , (6.10) 

where / and / are complex conjugate functions of z and z, respectively. By reparametriza- 
tion of the complex coordinates one may choose without lose of generality functions 
f(z) = z and f(z) — z so that their characteristic equation becomes A$ = z + z. A 
particular simple solution is, [25], 

e* = 7—^ • (6.11) 

{z + zf 

More general static solutions of type III are hard to construct beyond the ones that can 
be obtained from (6.11) by arbitrary reparametrization of the coordinates (z, z), whereas 
any dependence on t is dropped out for all practical purposes. 

Type D solutions are fully known and they include the Schwartzchild metric with 
mass parameter m as static example. In this class, there are additional relations 

d (e"*9(A$)) = , e"* (<9(A$)) 2 + 3md (e"*9(^$)) = , (6.12) 

and their complex conjugate, which are obtained from the condition 3\l/2^ / 4 — 2^3 2 
to different orders in the radial variable r. Integration of the first equation yields 
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exp(— $)<9(A$) = f(z;t), where f(z) is an arbitrary analytic function. Then, the 
Robinson- Trautman equation (6.4) simplifies to 

d (/e*) + 3md t e* = , (6.13) 

plus its complex conjugate relation. 

For / = the corresponding metric is t-independent and takes the form 

2r 2 / 2m \ 

ds 2 = ydzdz — 2dtdr — I K dt 2 (6.14) 

(l + fzz) V rj 

with constant K that can be normalized to or ±1. For K — 1 it coincides with the 
Schwartzschild metric in Eddington-Finkelstein frame, setting 

t — t s — r — 2mlog(r — 2m) (6.15) 

in terms of the usual time coordinate t s of the Schwartzschild solution. When / is not 
zero, one can perform coordinate transformation to make it constant, say / = —3m, in 
which case the function $ depends only on the combination z + z + t. Then, the geometry 
can be brought into the static C-metric form by further change of the coordinates (see, 
for instance, [28]). A notable example of this kind corresponds to the choice, [29], 

e* = ^— — , (6.16) 

(z + z + t) 3 / 2 v ' 

having 4mA 2 = 1. 

Type II solutions provide the last but most general class of space-times that are di- 
rectly related to time dependent solutions of the Calabi flow. They admit arbitrary initial 
data z; to) at some given time t = to, but unfortunately there are hardly any explicit 
solutions known to this day, apart from those obtained from type III configurations by 
allowing m/0, [25], as, for instance, the particular solution (6.11). Solution generating 
techniques have been subsequently employed toward the construction of more general 
type II space-times, but the results are rather inconclusive in the literature; see, how- 
ever, the recent work [30] on type II metrics which provides some explicit solutions that 
are not asymptotically flat. The main interest in type II solutions stems from the fact 
that both gravitational radiation and black hole formation can be studied simultaneously 
in exact terms in the vacuum. Thus, it is an outstanding problem in general relativity 
to explore the general structure of type II solutions; here come our results to make a 
decisive step in this direction. 



6.3 Further physical considerations 

It is clear from the discussion above that type III and type II Robinson- Trautman metrics 
provide the most general classes of gravitational metrics (6.1) in four dimensional space- 
times with parameter m zero and non-zero, respectively, whereas type N and type D 
metrics follow from them as special limiting cases. 
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Type II metrics provide a class of radiative space-times which tend to the Schwarzschild 
solution in the asymptotic future. Actually, approximate solutions of this kind were first 
examined in the literature on general relativity, [31], by perturbing the Schwarzschild 
metric so that the 2-spheres M defined at constant r and t represent distorted con- 
figurations in the multi-pole expansion of gravitational radiation. In particular, using 
the linearized approximation, as outlined in section 2.5, small perturbations associated 
to spherical harmonics with / > 2 were found to decay exponential fast by inducing 
quadrapole and higher order pole radiation, until the configuration eventually settles 
back to the round sphere as t — > oo. It was subsequently found that all dynamical type 
II solutions evolving toward the Schwarzschild solution exhibit Lyapunov stability using 
the quadratic curvature functional S(g), given by equation (2.16) for M = S 2 , [32]. As 
was also shown in section 2, S(g) varies monotonically with time and, in this case, the 
critical point where it assumes its minimal value is the round sphere that corresponds 
to the Schwarzschild space-time in Eddington-Filkenstein frame (6.14). In this case, 
the notion of extremal Kahler metrics that was introduced on general grounds by Calabi 
reduces to the constant curvature metric on S 2 that characterizes the end point of the dy- 
namical evolution. Finally, the (semi)-global existence and convergence of solutions with 
smooth initial data has been established in all generality, following arguments similar to 
the long time existence and convergence of solutions to the Ricci flow, thus establishing 
the universal limiting character of all asymptotically flat type II metrics, [12] (but see 
also [27, 33] for important earlier work on the subject). 

The global structure of Robinson- Trautman space-times has been carefully analyzed 
and extensions across the null hypersurface at t = oo have been considered in analogy 
to the Kruskal extension of the r > 2m Schwarzschild space-time. Two such solutions 
can be glued together along their Schwarzschild-like event horizon to form space-times 
that contain both a black and a white hole, as in the ordinary static case. However, the 
extension through the horizon is not smooth, in general, which implies that an observer 
living in this space can determine by local measurements whether or not he has crossed 
the horizon. These results were obtained by a detailed asymptotic analysis of the long 
time behavior of solutions to the Robinson- Trautman-Calabi equation, which upon trans- 
formation to Kruskal type coordinates were found to exhibit terms in the power series 
expansion with logarithmic coefficients, [34]. When the cosmological constant is not zero, 
the analysis of the solutions and their long time behavior essentially remains unaltered. 
However, the value of A enters in the space-time interpretation of the solutions and their 
long time asymptotic convergence to the Schwarzschild solution with a cosmological con- 
stant, e.g., the Schwarzschild-de Sitter metric when A > 0. An interesting feature of such 
generalization is that the continuation of the metric across the horizon can be made with 
a higher degree of smoothness when A > 0, [35], as compared to the more standard case 
with vanishing cosmological constant. 

Previous attempts in the literature to integrate the Robinson- Trautman equation for 
type II metrics include the application of prolongation methods, as in reference [36]. 
However, their results remain rather preliminary and inconclusive as they felt short of 
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expectation compared to other two-dimensional integrable systems that can be more 
easily treated by such tools. There are also some works on series solutions of type II 
metrics, [37], but a systematic understanding of their origin seems to be lacking. Our 
main contribution to this subject, according to the previous sections, is the ability to cast 
the class of all type II metrics into zero curvature form on the two-dimensional space (z, z) 
and construct systematically formal solutions based on free fields, via group theoretical 
methods. The only drawback of our construction is the inability to characterize the 
global structure of the resulting space-times, and their possible singularities, in terms 
of the free field data used in the Backlund transformations. Also, simple but explicit 
solutions of the equations are still unknown, as they may never exist simple models of 
spherical gravitational radiation in vacuum that only depend on finite number of moduli. 

Type III metrics, on the other hand, require separate investigation as they do not 
correspond to genuine dynamical solutions of the Calabi flow. True, they can be promoted 
to type II solutions of a very special type by allowing m 7^ 0, but they also remain 
largely out of reach apart from specific solutions, such as (6.11), and their descendants 
by coordinate transformations. The zero curvature formulation of type III metrics will 
be considered separately in the next section, where it is also explained how to extend 
the group theoretical methods of Toda theory in order to construct formal solutions in 
terms of free fields, as in the fully dynamical situation. We only note here that type III 
metrics are characterized by the non-vanishing component of the Weyl curvature tensor 
^3, which in turn implies that <9(A<3>) and <9(A$) differ from zero. Thus, type III metrics 
do not correspond to extremal Kahler metrics on S 2 , using the terminology introduced 
by Calabi, [6, 7]. On S 2 the extremal metrics are also constant curvature metrics with 
constant A$. 



7 The algebraic description of type III metrics 

A zero curvature formulation of the non-linear equation that characterizes type III 
Robinson- Trautman metrics, 

AA$ = (7.1) 

was given in reference [24] by taking suitable limit of the commutation relations of the 
algebra used for the Calabi flow. Here, we provide another zero curvature formulation 
of the same equation, which is not equivalent to that construction, but is more natural 
for applying the general framework of Toda field theories and their solution generating 
techniques as in the Calabi flow. 

First, notice that the equation (7.1) can be written in superfield form as 

ddT = ^ (7.2) 

using T = $ + 9^/ with bosonic components that do not not depend on time, i.e., t will 
be spectator if added by hand. Comparing the even and odd parts of the equation it 
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follows that 

d5$ = ^e* , 88V = , (7.3) 

which is equivalent to equation (7.1) after eliminating the field Thus, equation (7.1) 
admits a formal Toda field theory description with Cartan operator 8/88. 

Next, to appreciate the idea behind the new formulation of this equation that will be 
adopted in the sequel, let us briefly recall the essential features of Toda field equations 
based on the generalized system of commutation relations 

[X + (F), X~{Q)\ = H(S(FG)) , [H(f), H{Q)\ = , 

[H(F), X ± (G)] = ±X ± ((KT)Q) . (7.4) 

Lie algebras of this type provide a natural supercontinual extension of the more abstract 
class of continual Lie algebras (3.3) for any given pair of linear super-operators (K,S). 
Of course, as in the ordinary bosonic case considered in section 3, a supercontinual Lie 
algebra with characteristics (K, S) is isomorphic to that with (K = KS, 5 = 1), provided 
that the corresponding Cartan generators are related to each other by H(F) = H{SJ r ) 
when S is invertible. Thus, one typically makes the canonical choice (K, 1) without loss 
of generality. 

The associated system of Toda field equations follows from the zero curvature condi- 
tion [8 + A + , 8 + A-] = choosing the gauge connections 

A + = H(g) + X + (l) , A- = X~(e F ). (7.5) 

The resulting equations read, in general, 

8T = KQ , 8Q = Se* (7.6) 

and therefore one obtains 

ddT = KSe T = Ke T (7.7) 

after eliminating the field Q . Note that the final field equation only depends on the 
operator K and as such it is inert to the choice of algebra (K, S) or (K, 1). However, the 
intermediate equations (7.6), which transform into one another when different choices 
are being made, will be equivalent only if S is invertible operator. 

Thus, it does not make a difference whether one uses K — 1 and S = V T or the 
canonical choice K = T>t and S — 1, as we did before, for the zero curvature formulation 
of the Calabi flow. The situation changes drastically when the map H(!F) = H{SF) is 
degenerate, as in the case of nilpotent operators with S 2 = 0, for the two algebras with 
operators (K,S) and (K,l) are not anymore isomorphic. This particular situation is 
realized when S = 8/ 89 and the algebra with characteristics (8/88, 1) is not isomorphic 
to that with (1,8/88). Although the associated Toda field equation, which in this case 
reads as (7.2), does not know the difference, the intermediate equations (7.6) are actually 
sensitive on the choice. Hence, the natural question arises which is the most recommended 
choice of algebra in the present case. 
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Let us choose K — 1 and S = 8/89 and examine the detailed form of the equations 
(7.6), which in terms of components T = $ + 9^/ and Q = g + 9f become, respectively, 

= g , 8V = f , 

= ^e* , 9/ = . (7.8) 

Then, by eliminating / and g one arrives at the system (7.3) for the functions $ and \& 
that are required in all generality. On the other hand, the choice K = 8/89 and S — 1 
leads to a different system of intermediate equations (7.6), which in term of components 
read 

9$ = / , 8^ = , 

9 5 = e * , 9/ = ^e* . (7.9) 

Eliminating / and g, as before, we obtain the equations <9<9$ = ^e* and 8^ = 0, which 
also lead to AA$ = as required. However, the latter formulation is more restrictive 
than the first since d^f = 0, and its complex conjugate equation 9\& = 0, determine 
only a very special class of solutions to the second order equation <9<9$ = appearing 
in the system (7.3). Thus, generality requires making the choice (1,8/89) rather than 
(8/89, 1). The different choices yield equivalent descriptions only for constant curvature 
metrics in two dimensions, in which case the fourth order differential equation (7.1) 
becomes equivalent to Liouville equation. 

With these explanations in mind, we choose the super continual Lie algebra with 
characteristics K — 1 and S = 8/89, whose commutation relations take the following 
form in terms of components, using the decomposition of the generators A = A + 9Ai, 



[ffifoO, x ± (^)] = ±x±(^) , 

[h (<p), xtm = ±x±(^) , 

whereas the rest are trivial, 

[X+(v?), x -(^)] = o, 

[H (<p), X ± W] = , 



[X?(<p), X^m = ±H 1 (^) , 

[H^), Xfm = ±Xf(^) , (7.10) 

[x+(<p), xrm = o , 

[HM, H J (*P)] = 0. (7.11) 



These commutation relations should be contrasted with the choice made in reference [24], 
which is also quite general, but it constitutes an inequivalent algebraic description of the 
same problem which will not be utilized in the present work. 

The formulation that is adopted here is more appropriate for integrating the equa- 
tion AA$ = by group theoretical methods, at least formally, as in ordinary Toda 
systems. The present description has the advantage that t is completely decoupled from 
the commutation relations of the algebra, and, as a result, the smearing functions in the 
components of the generators, A ((p) and Ai(ip), can be completely dropped out. Put 
differently, the smearing superfields T can be considered as functions of 9 only, without 
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having any t dependence, so that they correspond to super-numbers tp + 9(pi. Also, the 
smeared generators can be simply defined by integration over the Grassmann variable 9, 

as 

A(F) = J A{9)F{9)d9 = Vl A + (7.12) 

and, as a result, the ip and ip dependence in the above commutation relations is super- 
fluous and factors out completely. 

Likewise, the highest weight state that is formally introduced to integrate the associ- 
ated Toda field equation (7.2) assumes a simpler form, compared to the Calabi flow, as 
it can be taken independent of t. Thus, here, we consider a normalized vacuum state 

\9 >= |0 > +9\0 >i (7.13) 

with defining relations, in smeared form, 

x+(f)\g >= o, <g\x~(F) = o, H{F)\g>=\pg> , (7.14) 

as before, where 

|;F>= J \9 > F(9)d9 = </?i|0 > +</?o|0 >i . (7.15) 
In terms of components these relations read 

#o|0> =0 , # |0>i= |0> , 

H^O > i= |0 >i for i = 0,1 , (7.16) 

and 

X +|0 >i= = X+\0 >i for i = 0, 1 (7.17) 

together with their Hermitian conjugate relations. 

Then, the group theoretical framework for Toda theories can be implemented in this 
case and the general solution of equation (7.2) takes the form 

F(z,z;9) =F (z,z;9) -log < 9\M-\z)M_(z)\9 > , (7.18) 

in accordance to the general expression (3.15) applied to supercontinual Toda systems 
with K — 1. JF = $ + ^^o, whose components are two-dimensional free fields 

%(z, z) = f(z) + f(z) , *o(*, z) = 4>(z) + $(z) . (7.19) 
Finally, the operators M± are given by the path-ordered exponentials 
M+(z) = Pexp (^J Z dz'e f{z,) [X+ + ij{z')X+]^ , 

M_(z) = VeKv(^j Z dz'e m [X~ +$(z')X-\ s j (7.20) 
that do not involve integration over the continuous time variable t. 
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Actually, there are several simplifications that occur in this case. Note that the 
commutation relations (7.10) and (7.11) contain a subalgebra generated by Xq 2 , Xf 
and Hi, which is sufficient to determine the structure of < Q\M+ 1 M-\Q >. The power 
series expansion of the path-ordered exponentials M± contain strings of these generators 
only, and, as a result, the evaluation of the corresponding vacuum expectation values, 
which is carried out as usual by pushing X + 's to the right and X~'s to the left, can 
never involve H and its action on the highest weight state \9 >. Since Hi\9 >= \9 >, 
as opposed to H \9 >= 9\0 > , it follows readily that the Grassmann variable 9 never 
appears as coefficient in the resulting expressions, and, for all practical purposes, the 
vacuum expectation value can be restricted to the state |0 > alone. Thus, 9 is also a 
superfluous variable in writing the general form of the solution, which simplifies to 

$(z,z) = %(z,z) - log ( < 0\M-\z)M_(z)\0 > ) . (7.21) 

At the same time, by comparing the odd terms of equation (7.18), we obtain to all orders 
the result 

y(z,z) = ^ (z,z) (7.22) 

thanks to the absence of ^-dependent terms in the expression o< OlAfj^M-lO >o- 

By construction, the power series expansion of a general field configuration $ in terms 
of free fields $o should coincide with the expansion that governs the Calabi flow when 
restricted to t-independent fields, since this limit is well defined within the general class 
of Toda systems. Thus, up to second order in powers of exp$ , it follows that for axially 
symmetric configurations with \l>o = a(z + z) + b and $ = c(z + z) + d, a broad class of 
solutions is 

e*° / 2a\ e 2 * / 2 4a T 7a 2 \ ^ ooX 

whereas ^ = ^ to all orders in the expansion, as required. In the special case of constant 
curvature metrics, which correspond to ^ = b, the system (7.3) reduces to the Liouville 
equation, 

dd§ = be* (7.24) 
and axially symmetric solutions of it assume the free field expansion 

* = *° + (Try" + (^ 2t " + "' < 7 - 25 ' 



according to the more general expression given above. 

Finally, it is interesting to identify the algebra whose local part is generated by X^ 1 , 
X^ 1 and Hi after dropping the superficial dependence on the smearing functions. Using 
the commutation relations (7.10) and (7.11), and changing base to 

60 = 7! (Xo+ + Xl+) ' 61 = 7! (Xo+ " Xl+) ' 

/o = i(I - + ir), fi = ^=(-X +X^) , (7.26) 
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it follows, setting also H 1 = h, that 

[e» fj] = hjh , [h, a] = d , [h, fi] = -fi . (7.27) 

These are precisely the defining relations of Kac's K 2 algebra, [38], which is the most 
elementary example of simple Lie algebra with infinite growth 4 . Equivalently, it can be 
thought as arising from the algebra 

[e% fj] = Sijhj , [hi, e 3 \ = K ijej , [h h fj] = -K^j) . (7.28) 

with i and j taking values or 1, using a degenerate 2x2 Cartan matrix with all its 
elements equal to 1. In such case, ho — h\ is a central element, as it commutes with all 
other generators e, and fi, and when it is factored out yields the K 2 algebra. 

The K 2 algebra was encountered in the literature before as the prolongation structure 
of the equation AA$ = that characterizes type III Robinson- Trautman metrics, [39]. 
It is rewarding that it also arises here using a somewhat different route of investigation. 
Our general construction based on Toda theories also connects it nicely with the algebra 
that determines formal solutions of the Calabi flow. It should also be noted that there is 
yet another way to write the type III metrics in zero curvature form, following [24] . The 
algebra used there is more complicated, as it depends on the additional parameter t - 
thought as spectator - in a non-trivial way, and it only contains K 2 in its zero modes. The 
intertwining of these two different algebraic descriptions of the same equation deserves 
further attention as it may provide the means to construct non-trivial conservation laws. 
This problem is left open for future work. 

8 Hierarchy of higher order flows 

Motivated by the general relation between the Ricci and Calabi flows, we introduce a 
hierarchy of higher order geometric flows by taking successive square roots of the time 
derivative operator. We will first describe the construction in two dimensions and outline 
an algorithm for casting the corresponding equations into zero curvature form. Finally, 
we will discuss the generalization of the hierarchy of geometric evolution equations to 
Kahler manifolds of arbitrary dimension. 

8.1 Higher order flows in two dimensions 

Let us assume that the conformal factor of the metric in two dimensions depends on 
infinite many time variables, U, as 

$ = *i, * 2 , *3, *4, •••) > 

4 The series of K n algebras is defined by the commutation relations (7.27) with i and j taking n 
different values. In all cases there is only one Cartan generator h and the growth of the algebras is 
infinite for n > 2. The case n = 1 corresponds to the si (2) algebra and naturally it is not included in 
this series. 
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and that there is an associated hierarchy of intrinsic geometric flows, 

w n = *"<*> • < 8 - 2) 

for appropriately chosen expressions of the scalar curvature R = — A$. The hierarchy is 
defined by imposing the following requirement for any two consecutive flows, 

<9e* <9 2 e* 

(8.3) 



dt n+ i dtl 



n 



This, in turn, implies d tn+1 & = — 9 t 2 n $ + (d tn <&) 2 , which suffices to determine re- 
cursively the form of curvature functional ty n (R), starting from the Ricci flow that 
corresponds to the choice = —R for the time variable t\. We have, in particular, 

*n + l = -^-afl, (8.4) 

which upon iteration yields \l/ n in terms of R. Also note that due to equation (8.3), the 
compatibility condition of the proposed hierarchy reads d tn+1 (\I/„exp$) = d tn (^ n +iexp^). 
Thus, in the appropriate sense, when the time derivatives act upon exp$ rather than <3>, 
the hierarchy (8.2) provides a system of independent flows; in this respect, it could be 
compared to other systems of non-linear differential equations with increasing order, such 
as the KdV hierarchy. 

This procedure can be simply implemented as stated above. However, there is an 
alternative description in two dimensions that will allow to prove the integrability of 
these flows by uncovering the relevant algebraic structures, as was successfully done for 
the Calabi flow. In particular, let us start from n-th flow and make the substitution 

W n -W n - e "aL' •-* + «.<*.*> (8-5) 

with the aid of a Grassmann variable 9 n , so that the (n + l)-flow is realized as super- 
evolution partner of the n-th flow. As a result, the order of the resulting parabolic 
equation is doubled, so that the n-th member in the family of flows (8.2) has order 2 n . 
The substitution made for $ can be understood as requirement of its super-analyticity 
on the super-coordinate T n = (t n ,9 n ) that extends t n in i? 1 ' 1 superspace; of course, 
9( n $ = ty n , which in turn determines ^ n +\ by working out the form of the equation in 
components. Also, consistency of the scheme requires that the Grassmann variables Q{ 
introduced in each step are all independent, i.e., OiOj + 6j6i = 25^. 

We have already seen in detail how this method yields the Calabi flow with respect 
to the time variable t-i (in the present notation). Further iteration yields the following 
eighth order differential equation, 

^ = A[(A$)(AA$) - AAA$] = A {BAR + AAi?) , (8.6) 
oh 
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using the expression ^ 2 — — AA$ for the Calabi flow. The integrability of this equation 
will be established later by implementing its algorithmic construction to the level of 
supercontinual Lie algebras. 

The structure of the equations becomes considerably more involved for higher values of 
n, as their order increases exponentially. We only give here the next one, for illustration, 

— = A[(A$) (A(^ 3 A$) - AAf 3 + 2^ 3 AA$) - (A# 3 )(AA$) 

ut^ 

- ^ 3 AAA$ + A (AAf 3 - A(^ 3 A$) - ^ 3 AA$)] , (8.7) 

where \l/ 3 = A((A<E>)(AA<E>) — AAA$). Higher order equations can be worked out at will, 
but their structure cannot be easily described in closed form. Unfortunately, we have no 
systematic way to iterate the procedure and obtain \l/ n for all values of n. Likewise, the 
choice of supercontinual Lie algebras that enables to cast them into zero curvature will 
only be prescribed algorithmically. 

Finally, one may formally append in the beginning of this list the heavenly equation, 
which is not parabolic but yields the Ricci flow by taking the square root of its time 
derivative operator. 



8.2 Integrability of higher order flows 

In order to understand how the hierarchy of higher order flows can be brought into zero 
curvature form, let us first consider the next example in the list given by the eighth 
order equation (8.6). In this case, the formal substitution d/dt 2 by d/d9 2 — 9 2 d/dt 3 
supplemented by the change of field variable $(z, z; t 2 ) to superfield J-{z, z; t 3 , 6 2 ) can be 
directly applied to the continual Lie algebra for the Calabi flow. Using its description in 
terms of components Hi((p) and Xf(ip) with % — 0, 1, as given in section 5, one obtains 
the following supercontinual Lie algebra with basic commutation relations 

[# (n x±(g)] = ±x±{rg) , [x±(n x^{G)\ = ±H {rg) , 

[#o(n xf(g)} = ±xt{rg) , ly+(.f), xr(Q)] = H^Q) , 

[h^F), xf(g)\ = T x±((v T r)g) , (8.8) 

whereas the rest are trivial, 

[H^r), x±(g)} = o , [x+(f), x-{g)\ = o , 

m?), H^g)] = . (8.9) 

Here, to simplify notation, we drop the indices from the new time variable t 3 and the 
associated Grassmann variable 9 2 and denote the corresponding super-evolution operator 
by V T , as before. 

Then, starting from the zero curvature formulation of the Calabi flow written in terms 
of components and promoting the smearing functions to superfields, so that the gauge 
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connections are taking values in the new super continual Lie algebra given above, we 
immediately arrive to the zero curvature formulation of equation (8.6). In particular, 
after some simplification of the original expressions used for the Calabi flow, we are led 
to propose the following choice of gauge connections, 

A + = # (&F) , 

A_ = X~(ddF) +X~(e :F ) , (8.10) 

where T and uj are now assumed to be i? 1 ' 1 mixed superfields with bosonic components 
that satisfy the relation 

V T u = -e- T ddT . (8.11) 

The zero curvature condition F z2 = [d + A + , d + A_\ = becomes equivalent to the 
following equation for the superfield J-{z, z; T) 

V T e T = -3d (e^ddF) , (8.12) 

which when expanded in terms of components T = $ + 9^f yields the system 

d t e* = dd (A* - ^A$) , = -AA$ . (8.13) 

Eliminating \1> leads to the eighth order evolution equation (8.6), as advertised. 

One can work out the form of the gauge connections, as well as the form of the 
supercontinual Lie algebra, in components by expanding everything in i? 1 ' 1 superspace. 
This will give rise to four bosonic components for each type of basic generators, which 
can be labeled as Hij(t) and X^(t) with % and j taking values or 1. The resulting 
expressions are quite lengthy and they are not included here. These generators, however, 
can be further promoted to elements of a supercontinual Lie algebra by changing d/dt to 
T>t and using i? 1 ' 1 superfields as new smearing functions. This procedure can be repeated 
indefinitely, and, as a result, one obtains a recursive construction of an infinite hierarchy 
of continual Lie algebras that are associated to the hierarchy of higher order flows. Thus, 
any member of the family is cast into zero curvature form by the appropriate choice of 
algebra. 

In general, a flow of order 2 n+1 , which is associated to a deformation variable t, is 
naturally obtained from a continual Lie algebra whose basic elements can be labeled as 
follows, 

H ili2 ... in (t) , Xt li2 ... in (t) , (8.14) 

with all indices taking values or 1. There are 2 n such independent elements for each 
type of generators, which follow by iterating the decomposition of i? 1 ' 1 superfields used at 
each level. Alternatively, they can be regarded as the bosonic components (in appropriate 
base) of three basic generators H(T) and X ± (T), which are labeled by a super continuous 
index in i? 1 '™ superspace with coordinates T = (t; 9i, 02, ■ ■ ■ , 6 n ), and likewise for their 
smearing functions. Thus, one may also consider the decomposition 

n 

A(T) = A (t) +"£8iMt) +J2°i8iMt) + ■■■ + ■ ■■0 n A 12 ... n (t) (8.15) 

i=l i<j 
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applied each basic generator H{T) or X ± (T). The coefficients in the expansion are fully 
anti-symmetric in their indices, and, therefore, the number of independent components 
for each generator also sums up to 2 n in this parametrization. 

It will be interesting to see whether the resulting hierarchy of infinite dimensional 
algebras can be put in more compact form in terms of H(T), X ± (T) so that the standard 
group theoretical methods of Toda field equations can be directly applicable to them for 
appropriate choice of Cartan operator. These steps will not be carried out explicitly here 
for the higher order equations but are left for the future. 



9 Generalization to arbitrary dimensions 

The hierarchy of geometric flows can be generalized to all Kahler manifolds of arbitrary 
dimension. Indeed, given the formal relation between the Ricci and Calabi flows, as 
outlined in all generality in section 2, we are led to propose a hierarchy of higher order 
geometric flows of the form 



imposing the requirement df n g ab = —d tn+1 g a i for all time variables t n . Everything that 
was said before generalizes quite naturally to all dimensions, apart from the underlying 
algebraic structures that are only useful for the zero curvature formulation of the flows 
on Kahler manifolds of complex dimension 1. We already have = —R ab and = 
d a dbR, and therefore all higher order equations assume the general form 

^-g al = d a d b V n (9.2) 

Ot n 

for appropriately chosen expressions V n (R). 

To illustrate the general structure of the equations, let us construct the next member 
of the hierarchy by computing the Calabi velocity of R, namely 

-^-R = -AAR - R al R al . (9.3) 

It is convenient at this point to introduce the fourth order elliptic operator of Lichnerowicz, 
which is defined for each Kahler manifold as 

dp = D LV - g ai R, a <p y - b . (9.4) 

According to the definition of Dl we have CR = AAR + R' ab R ab -, and, therefore, the 
next flow takes the form 

^-g al = d a B b (CR) . (9.5) 

It is a parabolic equation with respect to the new time variable t 3 and provides an eight 
order non-linear differential equation for the components of the Kahler metric. Similarly, 
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one may work out the form of higher order flows by iteration, but the final results are 
rather complicated. In general, the n-th flow is a parabolic equation of order 2 n with 
respect to the space coordinates. 

In all cases, apart from the (unnormalized) Ricci flow, the deformations are volume 
preserving, i.e., 

— ( dV(g) = I dV(g)AV n = (9.6) 

for n > 2. They also preserve the Kahler class of the metrics as for the Ricci and Calabi 
flows. In fact, for all these flows, the quadratic curvature functional S(g) that was intro- 
duced in section 2 exhibits an absolute minimum value equal to (J M R dV(g)) 2 / J M dV(g) 
that is attained by constant curvature metrics, when they exist on the Kahler manifold 
M. Also, the extremal metrics on M arise as local minima of S(g) for all higher order 
flows as well. 

More precisely, using the form of the flows d tn g a i = d a dbV n (R), we first compute the 
variation of S(g), 

9t n S(g) = ^- / R 2 [g]dV(g) = -2 / R(D L V n )dV(g) 

Ot n JM JM 

= -2(D L V n , R) = -2(V n , D L R) . (9.7) 

The result is obtained, as in section 2, by introducing the fourth order elliptic operator 
Di and removing the total derivative terms that arise in the course of the calculation. 
The Euler-Lagrange equation D L R = describes the extremal values of S(g) for all such 
flows, and the corresponding configurations are the extremal metrics on M, when they 
exist, satisfying the equivalent (but simpler) equation LR = 0, as for the Calabi flow. 
Actually, all these extrema are local minima because one finds that the Hessian form of 
S(g) in the directions t n and t m is 

= 2(D L D L V n , V m ) (9.8) 

OZ n Ob m 

when evaluated at the critical points. The operator Di commutes with its complex 
conjugate D L when a critical metric is being used and D L D L is a strongly elliptic eighth 
order operator that is self-adjoint and positive semi- definite. Then, for any given flow, 
the second variation is manifestly positive semi-definite, 

= 2(K, D L D L V n ) > , (9.9) 

implying that all extrema are local minima, as required. 

It is fair to say that the preceding analysis appeared for the first time in the original 
works of Calabi, [6, 7], who considered the general problem of deforming Kahler metrics 
on M within a given cohomology class, 

9ab = 9ab + d a d b U , (9.10) 
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using arbitrary t-dependent globally defined scalar functions u(z,z). The hierarchy we 
have presented above corresponds to specific choices of u, as given by the potential 
functions V n (z, z; t n ). The ordinary Calabi flow corresponds to the choice Vi{z, z; £2) = R 
and has the special property that S(g) decreases monotonically along it. All flows have 
critical points which are local (and possibly global) minima of S(g) associated to extremal 
Kahler metrics. We see that the only directions of deformation of a critical metric in which 
the expression (9.9) vanishes are those tangent to the orbit formed by all holomorphic 
transformations of M preserving the Kahler class of the metric. These transformations 
act as gauge group for the variational problem of minimizing the curvature functional S(g) 
for they leave it invariant. Conversely, any flow that deforms the metric transversely to 
the gauge group orbit has strictly positive definite second derivative, as given by equation 
(9.9). 

Finally, note that there are no curvature functionals in our disposal, at least for the 
moment being, which decrease monotonically along the corresponding higher order flows. 
It is an interesting variational problem to construct such expressions and study the long 
time behavior of their trajectories, as for the ordinary Calabi flow. Further implications 
of the higher order geometric deformations to problems of Kahlerian geometry are left 
open for future work together with their possible physical applications. 

10 Conclusions and discussion 

We have provided a unifying framework for formulating geometric evolution equations 
driven by the intrinsic curvature on Kahler manifolds, so that the metric deforms to 
canonical form within a given cohomology class. In two dimensions, in particular, these 
equations admit zero curvature description with gauge connections taking values in ap- 
propriately chosen infinite dimensional Lie algebras that incorporate the deformation 
variable t into their system. Such uneven treatment of the coordinates (z, z) and t al- 
lows to establish an algorithm for their formal integration in terms of a one-parameter 
family of two dimensional free fields, as in Toda field equations. At the same time, the 
evolution of any given initial metric exhibits dissipative behavior in t, since curvature 
perturbations around the canonical metric tend to decay exponentially fast. However, 
there is no contradiction of terms with the notion integrability in two dimensions, as t 
becomes internal index of the algebra of gauge transformations in the corresponding zero 
curvature conditions. Thus, it is possible to bring the equations into two-dimensional 
integrable form and hide all dissipative behavior in t in the internal space. 

The notion of supercontinual Lie algebras that was introduced here for the first time, 
has been particularly useful for establishing the zero curvature formulation of the two- 
dimensional Calabi flow by extending previously known results for the Ricci flow. This 
was implemented in practice by by extending t to i? 1 ' 1 superspace with coordinates (t, 9) 
and promoting the time evolution operator to super-evolution. Further iteration of this 
procedure gave rise to an infinite hierarchy of higher order geometric flows, which exist 
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on Kahler manifolds in all dimensions and they are integrable in two dimensions in 
the same sense that the Ricci and Calabi flows are. Actually, all these flows provide 
special examples of a more general variational problem posed by Calabi for deforming 
the metrics within a given Kahler class. In that context, one considers the quadratic 
curvature functional for a given Kahler manifold M of arbitrary dimension, whose local 
minima define the notion of extremal metrics on M, such as constant curvature metrics, 
provided that there are no obstructions for their existence. The general class of flows 
introduced by Calabi, d t g a i = d a dbU, exhibit common qualitative behavior in that they 
deform any initial geometric configuration toward the extremal metrics. The standard 
Calabi flow corresponds to the deformation driven by the scalar curvature, u — R, but 
there can be more arbitrary choices of the t-dependent scalar function u. The hierarchy 
we have constructed here selects some special functions u in a recursive way, such as R, 
£R, etc, and the corresponding higher order flows were found to admit zero curvature 
formulation in two dimensions. It will be interesting to examine whether more general 
Calabi deformations can admit similar description in two dimensions using appropriately 
chosen infinite dimensional Lie algebras depending on the form of u(z,z; t). 

The infinite dimensional structures that arise in the algebraic description of two- 
dimensional flows are quite novel, as they typically exhibit infinite growth. The prime 
example is the continual Lie algebra with Cartan operator K — d/dt used for the Ricci 
flow, which grows exponentially fast beyond its local part. The supercontinual algebra 
with K = T>t that was used for the Calabi flow is also known to have infinite growth, 
and the same behavior is expected from the infinite dimensional algebras associated to all 
higher order flows. Luckily, the zero curvature formulation of these flows involves gauge 
connections that take values in the local part of the algebra, and, hence, the complete 
structure of the algebra is not required for their Toda-like description and the formal 
construction of their general solution. However, it is natural to expect that the explicit 
construction of two-dimensional integrals requires better handle on the commutation 
relations of the algebras used in each case, beyond their local part, and the existence 
of invariant forms. It is conceivable that there will be no algebraic integrals of these 
equations, as in some other examples of exactly solvable systems whose zero curvature 
formulation is based on infinite dimensional gauge groups, like the long-studied Halphen 
equations, [40]. These issues remain out of reach at the moment, as more mathemati- 
cal work is required to bring the theory of such algebras, and their representations, to 
higher level of understanding. More generally, the very notion of integrability for two- 
dimensional systems based on infinite dimensional gauge algebras calls for further work, 
as there are only some isolated examples in the literature so far; see, also, references [41] 
and [42] where other examples are analyzed and compared to more traditional systems 
whose integrability properties are compatible with the Painleve criterion. 

The main physical application of the present work lies within the theory of spherical 
gravitational waves in vacuum in four space-time dimensions, as given by the general 
class of radiative Robinson- Trautman metrics. The formal integration of the Robinson- 
Trautman (two-dimensional Calabi) equation by group theoretical methods allows to 
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express all type II solutions (in Petrov's classification) using two-dimensional free fields, 
as in all Toda field equations. However, there are hardly any explicit solutions known to 
this day that may depend on a finite number of moduli of the deformed two-dimensional 
sphere sitting inside the four-dimensional radiative metrics. If that were the case, some 
simple solutions of gravitational radiation in the exterior of bounded sources would have 
been possible to manufacture and study in exact terms. Put differently, there seems to be 
no reasonable mini-superspace description of the Robinson- Trautman-Calabi equation, 
unlike the simpler case of two-dimensional Ricci flows that admit consistent truncations 
and lead to interesting axially symmetric deformations of the round sphere, as the sausage 
model. Mini-superspace truncations of non-linear evolution equations, when they exist, 
are also tractable by more direct computational methods that do not rely on the Toda field 
theory interpretation of the equations and the free field parametrization of their general 
solution. Thus, it appears that exact models of gravitational radiation depend on an 
infinite number of parameters and there is only a formal power series expansion for them 
based on the group theoretical interpretation of the general solution. The class of type 
III Robinson- Trautman metrics was also shown to define an integrable two-dimensional 
system with gauge connections taking values in an infinite dimensional algebra identified 
with Kac's K 2 algebra of infinite growth. 

Another interesting problem for future work is the generalization of the present for- 
malism to include geometric evolution equations on supermanifolds. So far, there has 
been no systematic work in this direction and it certainly deserves more attention. A 
seemingly related problem is the supersymmetric generalization of the Ricci and Calabi 
flows on Kahler manifolds. It will be interesting to examine whether in two dimensions, 
in particular, there are appropriate choices of continual and supercontinual Lie superal- 
gebras that can accommodate such generalizations in the form of zero super-curvature 
conditions. A different but illustrative example is provided by the N — 1 supersymmet- 
ric generalization of the heavenly equation, whose bosonic version is viewed as continual 
Toda system in two dimensions with gauge connections taking values in the local part of 
the SU (oo) algebra. In the supersymmetric extension of the heavenly equation, [43], one 
considers the Lie superalgebra sl(N\N + 1) in the large N limit in which the system of 
simple roots is odd and corresponds to the super-principal embedding of osp(l\2) in it. As 
a result, the Cartan operator of the corresponding continual superalgebra is d/dt, equal 
to the square root of the Cartan operator of SU (oo), and it happens to coincide with 
the Cartan operator of the continual Lie algebra used for the Ricci flow. This occurrence 
might be only the beginning of a systematic pattern for constructing supersymmetric 
extensions of the hierarchy of higher order geometric flows, together with their zero cur- 
vature formulation, by also turning ordinary derivatives into supersymmetric covariant 
derivatives in two dimensions. It remains to be seen how far this construction can be 
pushed and what is the systematics behind it. 

The use of infinite dimensional algebras that incorporate the time variable t into their 
defining system of commutation relations might be of more general value for the reformu- 
lation of many dynamical problems other than the geometric evolution equations. If this 
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possibility materializes it will certainly provide new ways to address difficult dynamical 
problems in terms of new classes of infinite dimensional Lie algebras and their represen- 
tations. It may also fertilize the interplay between physics and mathematics in an area 
where Lie algebras of infinite growth become relevant and demand better understanding 
of the whole subject. 

Finally, it is quite interesting to consider the interplay between physics and mathe- 
matics for other classes of geometric evolution equations and develop new tools for their 
study. In this context, a unified picture of intrinsic as well as extrinsic flows could emerge. 
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